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Abstract 

Low-energy results from measurements of leptonic dipole moments are used to derive 
constraints on the CP-violating phases of the dimensionful parameters of the minimal 
supersymmetric extension of the standard model (MSSM). We use these (known) bounds 
to investigate the impact of these phases on CP-even cross sections at high-energy e^e~ 
and e~e~ colliders. To that end we define two measures of the significance with which the 
existence of non-vanishing phases could be deduced from the measurements of these cross 
sections. We find that highly significant evidence for deviations from the CP-conserving 
MSSM could be obtained at the next e~^e~ collider even if the electric dipole moment 
of the electron is very small or zero. We also analyze a CP-odd final state polarization, 
which can be large when two different charginos or neutralinos are produced. Finally, we 
study correlations between the phase-sensitive observables. 
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1 Introduction 



CP violation was observed first in the neutral kaon system [1], and has recently been found 
in 5— meson decays [2]. In addition, CP violation constitutes one of the conditions for a 
dynamical generation of the cosmological baryon asymmetry [3]. In the Standard Model (SM), 
which contains only one physical neutral Higgs boson and assumes neutrinos to be massless, 
the only source of CP violation is the complex phase of the quark mixing matrix [4] .* 

Supersymmetry (SUSY) is now widely regarded to be the most plausible extension of the 
SM; among other things, it stabilizes the gauge hierarchy [6] and allows the grand unification 
of all known gauge interactions [7]. Of course, supersymmetry most be (softly) broken to be 
phenomenologically viable. In general this introduces a large number of unknown parameters, 
many of which can be complex [8]. In the most general minimal supersymmetric standard 
model (MSSM) 44 phases cannot be removed by suitable redefinitions of fields and remain as 
"physical" phases in the model. For example, they have a direct impact on the mass spectra 
as they enter most mass matrices in the Lagrangian. Of course, one can use more specific 
assumptions on the soft breaking terms and/or an underlying GUT theory to get simpler 
versions of the MSSM with a smaller number of parameters, but the price for doing so is the 
loss of generality. 

CP- violating phases associated with sfermions of the first and, to a lesser extent, second 
generation, and with the chargino/neutralino sector, are severely constrained by bounds on 
the electric dipole moments of the electron, neutron and muon. However, as emphasized in 
[9, 10, 11, 12] cancellations between different diagrams allow some combinations of these phases 
to be rather large even for a sfermion mass spectrum accessible at the expected center-of-mass 
energy of a possible next linear e^e~ collider (LC). Even in models with universal boundary 
conditions for soft breaking mass at some very high energy scale, the relative phase between the 
supersymmetric Higgsino mass parameter /i and the universal trilinear soft breaking parameter 
^0 can be 0{1) [13]. If universality is not assumed, the relative phase between the U{1)y and 
SU (2) gaugino masses may also be large. 

In the past few years a lot of effort has been devoted to analyses of the physics output 
that can be expected from experiments at the LC, including a possible e~e~ option [14]. Work 
towards the design of such a device has also made great progress. Today it is assumed that 
it will (initially) have a center-of-mass energy ^/s in the range between 500 GeV and 1 TeV, 
an integrated luminosity of at least several hundred fb~^, and adjustable polarization for both 
beams. Detailed analyses [14, 15] have established that sparticles with mass can easily 

be discovered at an LC. Moreover, many of their properties (masses, spins, some couplings) 
can be measured precisely. 

Unfortunately most of these analyses [15] show the dangerous tendency to neglect phases, 
which are actually free parameters of the model and are not necessarily neghgibly small. Note 

*The observation [5] of neutrino flavor oscillations opens the possibility that the ncntrino mass matrix 
contains non trivial CP violating phases, but this has not yet been confirmed experimc^ntally. In principle CP 
could also be violated in the SM by the QCD O—term, but bounds on the electric dipole moment of the neutron 
force |6Iqcd| to be < 10''^°. 
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that both masses and couphngs depend on these phases, which will hence have a direct impact 
on sparticle production cross sections and decays. Neglecting non- vanishing phases when deter- 
mining real parameters from experimental data could thus lead to wrong inputs for attempts to 
reconstruct the underlying theory at the unification scale. On the other hand, the construction 
of sizable and experimentally accessible CP-violating observables is rather difficult in most of 
the production channels at e"'"e~ colliders, as at least one secondary decay has to be included 
in the analysis. At the tree level nonzero CP-odd asymmetries can only result if the decaying 
particle has nonzero spin, which should be at least partly reconstructed from its decay products. 

We therefore first perform a rather general analysis of the impact of non-vanishing CP-odd 
phases on CP-even cross sections. We work in the framework of the MSSM with non-vanishing 
CP phases. We assume flavor universality for soft breaking terms associated with sfermions of 
the first and second generation, but we do not assume any specific model for SUSY breaking. 
Our free parameters are specified at the typical energy scale of an LC. The basic idea of this 
work is to take today's low energy data, such as lower mass bounds and bounds on leptonic 
dipole moments (rfg and a^), as constraints for a parameter space scan. We then use the 
resulting, low energy compatible points to check whether high energy experiments at an LC (in 
either the e'^e" or the e~e~ mode) could provide additional information on phases. We restrict 
ourselves to the following total, unpolarized cross sections: 

i,j = l,...,4; (1.1a) 

t,j = 1,2; (1.1b) 

t,j = l,2; (1.1c) 

i,j = 1,2. (l.ld) 

There is a complementarity between the leptonic dipole operators and the high energy pro- 
duction amplitudes. Since several diagrams involving neutralinos as well as charginos contribute 
coherently to the low-cncrgy observables, they can only give bounds on combinations of phases. 
In contrast, high energy observables can be used to investigate the different sectors of the theory 
separately. As our aim is to study the impact of low energy compatible, non-vanishing phases 
on the cross sections, we assign a significance S{f 1/2) to each final state, defined as difference 
in production rates between a CP-conserving point (CPC-point: real parameters, all phases 
identical to zero or vr) in parameter space and a CP-violating point (CPV-point: same absolute 
values of parameters, but non-vanishing phases) normalized to the statistical error of the cross 
section in the CPC-point. Since the phase dependence of a given cross section might partly 
arise from kinematical effects (kinematical masses depend on the phases), we also introduce a 
second significance 5(/i/2), where the CPV-point is chosen such that the masses of two neu- 
tralinos and one chargino coincide with the CPC-point; this can be achieved by adjusting the 
absolute values of the relevant dimensionful input parameters. 

We find that these significances can be very large for some reactions of the types (1.1a) and 
(l.ld), but are usually small for (1.1b) once the low-energy constraints have been taken into 
account. Moreover, if the absolute values of the input parameters, or three chargino/neutralino 
masses, are kept fixed while the phases are varied randomly, there is no visible correlation 
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between these high-energy significances and de- On the other hand, in most cases sizable high- 
energy significances are strongly correlated with each other, and slightly less strongly correlated 
with a^. 

Strictly speaking these significances only measure deviations from the CP-conserving version 
of the MSSM. These deviations might also be explained by some extension of the MSSM without 
invoking new sources of CP violation. If some deviation from the CP-conserving MSSM is 
observed, a more direct probe for CP violation in the production and decay of superparticles 
thus becomes important. We therefore compute the CP-odd polarization of charginos and 
neutralinos that is normal to the production plane, and find that in all cases considered, it 
can be sizable for neutralinos; a large CP-odd polarization of charginos is possible only at 
large and small tan/3. In contrast to earlier, related work [16, 17] we emphasize a detailed 
semi-analytical understanding of the observed effects; isolate the measurements that hold the 
most promise; and analyze the correlations between various phase sensitive observables. 

The remainder of this paper is organized as follows. In Section 2 we briefiy review the 
mass spectra and mixing patterns of the sleptons, charginos and neutralinos. This section 
also contains an overview of the relevant parameters. After summarizing the relevant parts 
of the MSSM Lagrangian in Section 3, we present in Section 4 the analytical expressions for 
the SUSY contributions to de and a^^, and discuss briefly possible scenarios for suppressing 
these leptonic dipole moments while keeping some phases sizable. This Section also discusses 
numerical constraints on these phases in three benchmark scenarios where selectrons as well as 
the lighter neutralino and chargino eigenstates can be produced at a 500 GeV e"'"e~ collider. 
Section 5 summarizes the well-known results for total cross sections of the production channels 
(1.1). We also give results for the components of polarization vectors for reactions (l.la,b). 
The significances are introduced in Section 6. In Section 7 we show the most important results 
of our detailed numerical analysis of the high energy observables. Section 8 completes our work 
with a brief summary of our findings and some conclusions. 

2 Particle mixing 
2.1 Slepton mixing 

As mentioned in the Introduction, we will assume that flavor mixing is negligible in the slepton 
sector. This can e.g. be motivated by the very tight experimental constraints on branching 
ratios for lepton flavor violating decays like n — > 67, — > 3e etc. The simplest way to satisfy 
these bounds on flavor changing processes is to assume that soft SUSY breaking parameters 
in the slepton sector arc the same for the first and second generation, as is the case in most 
models that attempt to describe SUSY breaking by a small number of parameters (which are 
usually defined at a high energy scale). The only relevant mixing in the slepton sector then 
occurs between SU{2) doublet sleptons II and singlets Ir. The squared mass matrix A4j in the 
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basis {II, Ir) is given by [18] 



A^H ( y! ) . (2.1) 



The elements of this matrix are defined as 



Xi = + m?^ + ^ (M| - 2M^) cos 2/5, (2.2a) 

Yi = + m?^ + (M^ - M|) cos 2/?, (2.2b) 

\Zi\ = m^lAf + //tan/3|, (2.2c) 

arg(Zf ) = 0^- = arg (- A - tan (3) , (2.2d) 

where m; is the mass of the charged lepton m?^ and Ai are soft SUSY breaking parameters, 
which we assume to be the same for the first and second generation, /i is the Higgsino mass 
parameter, and tan /3 is the ratio of vacuum expectation values (vevs) of the two neutral Higgs 
fields. In general, /i = |/i|e*''^'" and Ai = \Ai\e^'f'^ can be complex, while all other parameters 
appearing in eqs.(2.2) are real. 

A4j can be diagonalized by a unitary transformation 

C/lA^?C/^ = diag(m?,m|), (2.3) 

with the mass ordering < m? by convention. The diagonalization matrix Uf can be param- 
eterized as 

f cosOi -sin%e-^^f\ 
' \sin%e^^f cos^;- J' ^ > 

where — 7r/2 < 9i< 7r/2 and < 0;- < 27r. Defining 



Mf ^ ^^—^ = (2.5a) 



A,- ^ ml - ml = sJ{Xj- Y^^ + A\Zf, (2.5b) 
the slepton mass eigenvalues and mixing angles are given as 

sin 26, = -2^; cos 2^.- - ~ "^^ (2.6b) 
' ' A,- 

Eqs.(2.6b) and (2.2c) show that slepton left-right mixing is suppressed by the corresponding 
lepton mass, but is enhanced for large tan/3 and large \^\. 

As sncutrinos are only present as components of left handed superfields in the MSSM, there 
is no partner to mix with and the mass simply reads 

m| =m?^ + ^cos2/3M|. (2.7) 
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2.2 Chargino mixing 

The Dirac mass matrix for charginos mixes the SU (2) gaugino and the charged Higgsinos 
h"^ . In the basis (w^, h~) it is given by [18] 

•^^^^ V2M,^sin/3 „ J' ^2-^) 

where the soft breaking mass parameter M2 for SU{2) gauginos is taken to be real and positive; 
this can be achieved without loss of generality by appropriate field redefinitions. This complex 
mass matrix is asymmetric and hence has to be diagonalized by a biunitary transformation 

UrMcUI = diag (m-±, m-±) , (2.9) 

with the mass ordering m^± < m^± as convention. The mixing matrices may be written as [19] 



Ul={ ''"'^"f' ), (2.10a) 

' —smcpLe^P^ cos (pL ' 



with — 7r/2 < (t)L,R < 7r/2 and < 71,2, Pl,r < 27r. 71 and 72 denote two possible Dirac phases 
which have to be introduced to ensure that the mass eigenvalues oi M.c ai'e positive and real. 
The parameters of Ul and Ur can be determined from A^^^A^c and M-cMc, respectively. 
Introducing the quantity 



Ac = { (M| - + 4M^ cos'(2/3) + 4M^ (M^ + |/x[2) + 8M^[/x| cos 0^M2 sin 2/?} 



1/2 



= (2.11) 
the squared mass eigenvalues are: 

mL =^(M2+|;x|2 + 2M^tAc), (2.12) 

while the mixing angles can be computed from 

-M,^+ 1/^1' + 2M^ cos 2/3 
cos20i = , (2.13a) 



A, 



c 



sin20i = ^^^^ {M^cos'^p+\ii\^sm^P + M2\ii\cos(j)^sm2(3Y^^ , (2.13b) 

-M| + |//[^-2M^cos2/3 

cos20i^ = , (2.13c) 

sm2(j)R^ -2A/2Mty (Ms^sin^/? + [/xl^cos^/? + Msl/^l cos0^sin2/5)^^^ (2.13d) 



and the phases are 



= I T'T r«M - P.14a) 

1/^1 COS 0^ + cot pM2 

\a\ sin (hn . s 

tan f3R = ^ ; (2.14b) 



M^|//|cos0^sin2/3 + M2 (^m^± - |At|^) 



I // 1 sm 0^ sm 2p 



cos0^ (^m^i - Af^j + M^M2sin2/3 



cot72 = ^ ^ . (2.l4d) 



I// 1 sin0^ (^m^± - M|) 



2.3 Neutralino mixing 

The neutrahno mass matrix A4n mixes the neutral components of both Higgsinos h^^^ with 
hypercharge ±1/2, the f^(l)y gaugino B and the neutral SU{2) gaugino W3. The mass matrix 
in the basis {B, W3, h% hi) reads [18] 

/Ml 



Mn = 



— Mz cos /3 sin sin /3 sin \ 

M2 M^cos/Scos^w — M^sin/^cos^w 

— cos /3 sin Mz cos /9 cos 6*-^ 

y sin /5 sin 6'w — sin/3cos^w 



-// 
-/X / 



(2.15) 

The U{1)y gaugino mass parameter Mi = |Mi|e"^i is in general complex. This symmetric mass 
matrix is diagonalized by a unitary transformation 

N'^Ai^N = diag(m^o, m^o, m^o, m^o), (2-16) 

i.e. the n— th mass eigenstate* is given by the complex conjugate of the n— th column of A^. 
Although the neutralino mass matrix can be diagonalized analytically even for complex param- 
eters [20], the general expressions are too lengthy to reproduce here. Of course, a numerical 
computation of N is straightforward. However, in order to qualitatively understand mixing in 
the neutralino sector, a perturbative diagonalization of the mass matrix (2.15) is often suffi- 
cient. Here Mz is considered to be a small parameter compared to |Mi|, M2 and \fi\. Keeping 
all terms up to first order in Mz, as well as a few C(M|) terms that will be important later, 
one finds for the masses and eigenvectors: 

m^o ~ |Mi| + Srm, X? ^ e^^^/^ ^^^^ ^^^^ ^^^^ ^2.17a) 
m^o ~ M2 + 5m2, ~ {521, 1, S23, ^24) /N2; (2.17b) 



*Whcn written as a row vector in the {B, W3, /i", /t") basis, the mass eigenstate satisfies AdNiXnV 
m^o {Xn)^ ('^o sum over n), i.e. it is not an eigenvector of Mn in the usual sense. 
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i(</.^+7r)/2 



^x°3 - l/^l' xi - — ^ — (<^3i, fe, 1, 1) ; (2.17c) 

- l/^l' X4 - (<^4i, S42, 1, -1) . (2.17d) 

In eqs.(2.17) we have assumed the ordering \Mi\ < M2 < If these three mass parameters 
are ordered differently, the eigenstates in eqs.(2.17) are no longer labeled in order of increasing 
mass. Note that the eigenvalues do not receive 0{Mz) corrections. However, mixing between 
gauginos and Higgsinos is generated at this order, and B — W3 mixing is generated at order 
M|. These mixings are described by the complex quantities Sij in eqs.(2.17); they are given 
by: 



<^i2 = 7172 rTTTKTn^. HTTH^ ' (2.i8aj 



M| sin Ow cos Ow [\Mi\ 




{Mi- 




2 _ 


|Mi|2) 



Mz sin OwiM* COS P + n sin (3) m -to-u\ 

~ |/iP-|MiP ' ^^-^^^^ 

. Mzsin^w(M*sin/3 + /xcos/5) ,^ „ , 

H2-|Mi|2 ' (2-^^") 

M| sin 9w cos 9w [M| + M1M2 + sin(2/5) {M2I1* + Mi/^)] 



M2Cos6'w (M2COS/5 + /xsin/5) /o 10 ^ 



Mzcos6'w (M2sin/5 + /icos/3) /oiofN 

'^^^ H^^Mf ' ^ ^ ^ 

M^sin^wlsin/i- cos/3) /oiqn 

'^ai = (|^|2_|M,|2) (^1 - H ; (2-18g) 

= ^^^osMsin/3 - cos/3) _ . 
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. _ Mzsingw(sin/3 + cos/3) 

(^41 (|^|2_ |Mi|2) (VWi + ) , (2.181) 

MzCOsew(sin/3 + cos/3) 1Q•^ 
^^42 = (|^|2_M|) + ■ ^2.18j) 

A?"!, in eqs.(2.17a,b) are normahzation constants, which differ from unity at 0(M|). Note 
that the phases of the 0— th order eigenstates have been factored out in eqs.(2.17); this gives 
more symmetric looking expressions for the 5ij. Finally, the 0{Mz) mass shifts of the gaugino- 
like states are given by 

'^^i = - | if u7i2 Il^il + l/^l sin2/3cos(0i + 0^)] ; (2.19a) 



- - I , J [M2 + Sin 2/? COS 0^] . (2.19b) 



|/i|2 - M; 
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Eqs.(2.18) and (2.19) show that the expansion wiU break down if — |Mi| or — M2 
becomes close to Mz in absolute value. In other words, even if the unknown dimensionful 
parameters in the mass matrix (2.15) are all » Mz, there can still be strong Higgsino-gaugino 
mixing if some of these parameters have similar absolute values. 

2.4 Relevant parameters 

The mixing patterns in the part of the SUSY spectrum which will be relevant for the remainder 
of our work depend on 10 SUSY parameters (plus some SM parameters whose values are already 
known accurately). Some of these parameters {rni^. fni^, \Ai\, (f)A, \Mi\, 0i) only enter in a 
single sector (sleptons and neutralinos, respectively), while M2 appears in both the chargino and 
neutralino mass matrices, and 0^ and tan/3 affect all three sectors. Therefore the mixing 
patterns in the separate sectors are partly correlated to each other. In particular, choosing 
the parameters of the neutralino mass matrix completely determines the chargino mass matrix 
as well. Moreover, increasing suppresses gaugino-Higgsino mixing, but enhances II ~ Ir 
mixing. Finally, taking tan/3 ^ 1 again enhances Il — Ir mixing, but reduces the impact of all 
phases on the physical masses. 

3 Interaction Lagrangian 

In order to make our paper self-contained, and to fix the notation, this section is devoted to 
a short collection of the relevant pieces of the interaction Lagrangian expressed in terms of 
physical mass eigenstates. 

3.1 Interactions involving SM gauge bosons 

First of all, the well-known SM coupling between charged leptons and gauge bosons is given 

by 

^ih,z = elY {A^Q^'^Pa + Z^QfP^) /, (3.1) 

where e is the QED coupling constant, and Pa, ol £ {+, — } = {-R, L}, are standard chirality 
projection operators, defined as 

1 ± 75 , . 

P± = (3.2) 

The linear charges in eq.(3.1) are 

g+'' = = 1; (3.3a) 



Qz = —n r- - ; (3.3b) 

sm^wcos^w V 2 ' 
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tan 



(3.3c) 



Slepton L ^ R mixing does not affect the couplings between sleptons and photons. Moreover, 
in case of selectrons L — R mixing can safely be neglected for high-energy applications. The 
vertices with two charged sleptons and one gauge boson are defined via the momentum-space 
Lagrangian 

= e [A,Qf + Z,Qf) {h + k,Y~k{h)l,{k,)\ (3.4) 



where i,j e {i?, L}. The corresponding linear charges Q^f'z are 



— ^ij ] 



1-3 



tan Q 



1 



w 



2 cos sin 

The couplings between physical Majorana neutralinos and the Z boson are given by 



2 cos sin 



where the linear charges Q^o are defined as 



(3.5a) 
(3.5b) 

(3.6) 
(3.7) 



The first equality in eq.(3.7) follows from the Majorana nature of the neutralinos. Of course 
there is no neutralino-photon coupling. 

Finally, the interactions between neutral gauge bosons and charginos are given by 



with 



-1 



(sin^ e^^Sij - {W±)ij) . 



cos sin 

The matrices [W±)ij can be obtained from the chargino mixing matrices via 

and read explicitly in terms of chargino mixing angles and phases as 

3 



4 + icos20i -i sin 201.6-*^^ 
-isin20ie^^^ |-icos20L 



•isin20Re"^(Ti-''«-T2) 



-isin20i?e*(^i-'^«-^2) 



cos 2(f) 



R 



(3.8) 

(3.9a) 
(3.9b) 

(3.10) 

(3.11a) 
(3.11b) 
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3.2 Slepton interactions with a chargino or neutralino 

The neutralino-slepton-lepton vertices receive contributions from both gauge and Yukawa in- 
teractions: 

^U^o = ' r(Gg + Y^^) P^x% + h.c, (3.12) 
V2sin^w 



with 



Gr. ^ -2imie^^Nu{Ui)2j; (3.13a) 

G± = (tan^wiVi* + ^i) {UiW, (3.13b) 

= -^/2YlN,,{U^)^J■, (3.13c) 

Y^ = -v^l^7V3*.(C/-)2,.. (3.13d) 



Here the dimensionless, rescaled Yukawa couphng Yi is given by 

Y, = . (3.14) 

Note that we have to keep terms oc Yj, as well as a non-trivial sleptonic mixing matrix Ui, when 
computing leptonic dipole moments. On the other hand, for high-energy applications Y^ can 
be set to zero, which implies Y-"^ = in case of selectrons. In the same limit L — R mixing can 



be neglected, in which case the gauge contributions simplify to 

G~j = — 2 tan O^fiiNiiSRj; (3.15a) 
G± = (tan^wA^i* + Slj. (3.15b) 

The couplings between chargino, sneutrino and lepton also receive gauge and Yukawa con- 
tributions: 

where 

Na,i = -SaL{UR)a + SMULhYi. (3.17) 

As before, the term oc Yg in eq.(3.17) can be dropped in collider physics applications, but it 
has to be kept when computing the leptonic dipole moments. 



4 Low energy constraints 
4.1 Experimental constraints 

In this paper we are only interested in purely leptonic processes. We therefore ignore the 

(quite stringent) experimental constraints on the electric dipole moments of the neutron and 
mercury atom. The main reason for this choice is that leptonic processes suffer much less 
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from uncertainties due to non-perturbative strong interactions. For example, ref.[21] finds that 
different models relating electric dipole moments of quarks to that of the neutron or Hg atom 
differ by typically a factor of two. Since large phases in the hadronic sector can be tolerated 
if there are cancellations between different contributions, which have different hadronic matrix 
elements, a conservative interpretation of the experimental bounds on dn tends to give [11] 
significantly weaker constraints on model parameters than the bound on the electric dipole 
moment of the electron does, even if one assumes some connection between the CP-violating 
phases in the squark and slepton sectors. The only CP-violating (more exactly, T-violating) 
low-energy quantity of relevance to us is therefore the electric dipole moment of the electron 
dg.. Given our assumption of flavor universality of the soft breaking terms in the slepton sector, 
at least as far as the first and second generation are concerned, the bound on the electric 
dipole moment of the muon [22] need not be considered separately: since {di)svsY oc nii, all 
combinations of parameters that satisfy the constraint on the SUSY contribution to de will be 
at least five orders of magnitude below the maximal allowed SUSY contribution to rf^. 

On the other hand, our assumption of universal sleptonic soft breaking terms for the first 
two generations also implies [23] that the measurement [24] of the anomalous magnetic moment 
of the muon, = (^f^ — 2)/2, gives a ^i^f/iter constraint on SUSY parameters than the anomalous 
magnetic moment of the electron does. The reason is that for universal soft breaking masses 
the SUSY contribution to these leptonic magnetic moments is essentially proportional to the 
squared mass of the lepton, and the experimental errors satisfy [22, 24] 5a^/m^ < 6ag/ml. The 
second low-energy quantity of relevance to us is therefore a^. Note that the SUSY contributions 
to ttfj, and de show very similar dependences on the absolute values of the relevant parameters; 
however, de receives nonvanishing contributions only in the presence of nontrivial phases, while 
the contribution to |a^| becomes maximal if all phases are zero or tt. 

The SM prediction for dg is negligible. The current experimental measurement [22] 

((^e)exp = (0.069 ± 0.074) X 10-2^ e • cm (4.1) 

can therefore directly be translated into a 2a range for the supersymmetric contribution to de- 

- 0.079 X 10"^*^ e • cm < (4)susy < 0.217 x 10"^*^ e • cm. (4.2) 

The interpretation of the most recent measurement [24] of a^, 

(aM)exp = (11659208 ± 6) x 10"^°, (4.3) 

is less clear. The reason is that non-perturbative hadronic terms do contribute to a^^, at about 
the 10~^ level. In principle this contribution can be calculated from experimental data using 
dispersion relations [25, 26]. Unfortunately calculations based on different data do not quite 
agree, although the discrepancy has become smaller after the recent release of corrected data 
by the CMD-2 collaboration [27]. Using e+e^ annihilation data as input tends to give an SM 
prediction which falls a little short of the experimental value (4.3). A recent analysis which 
includes all existing e+e" data [28] finds 

(a,.)sM = (11659180.9 ± 8.0) x 10"^°. (4.4) 



12 



Adding all errors in quadrature, this gives a ~ 2.1a discrepancy. On the other hand, using data 
from r decays gives [28] 

K)sM = (11659195.6 ± 6.8) x 10"^°, (4.5) 

which is only ~ 1.4cr below the measurement (4.3). Since even the e^e~ data lead to a less than 
3cr discrepancy between the prediction for and measurement of a^, we do not want to claim 
evidence for a non-vanishing SUSY contribution. In order to be conservative, we construct 
the upper limit of the "2(T allowed" range for (a^)susY = ('^/^)cxp — ('^/x)sm by using the lower 
value (4.4), reduced by two combined standard deviations, as our estimate of (a^)sM- Similarly, 
the lower end of this "2(T range" is obtained when (a^)sM is estimated by adding two standard 
deviations to the higher value (4.5). This gives: 

- 5.7 X 10-^° < (a^)susY < 47.1 x 10"^°. (4.6) 

The upper bound in (4.6) constrains the SUSY parameter space only for large values of tan /3, 
but the lower bound is significant also for moderate tan (3. 

4.2 Analytical results 



'7 



/7 i~ 




Figure 1: SUSY contributions to leptonic dipole operators 

The supersymmetric one-loop contributions to lepton dipole moments are shown in Fig. 1. 
The left diagram depicts the neutralino contribution while the right one contains the chargino 
contribution. Using the interaction Lagrangians given in Section 3 we find for the chargino 
contribution to the electric dipole moment of the electron 

The chargino loop contribution to the magnetic dipole moment of the muon is 

,2 



2 ( 

zb X I ^Tfi Tfi^ 

«m)susy = E \ :^fii^imcuc%) + -JLf,{x,){\cu\' + \cm\') \ • (4i 

i=l I Xi xf 
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The corresponding results for the neutrahno contribution read 



-1 



4 2 



e I SUSY 967r2 --^ m^o 



(4.9) 



(«m)susy = X X ] ^/2(l/ia)3f^e(n^i„nRi„) + "^hiViaWnuc^ + 



The variables Xi and i/j^ are defined as 



2 

2 ' 



2 

m~ 

Via 2~' 



(4.10) 



(4.11) 



and the loop functions fi are 

m = 2{z^-l)^ ^'' -Az + 3 + 2\ogz); 

h{z) = 7r^^(^^ - + 3z + 2 + 6^1og^) . 



(4.12a) 
(4.12b) 
(4.12c) 



These functions arc normalized such that /j(l) = 1, i = 1, 2, 3. Finally the coupling coefficients 
cai and nAia can be written as: 



Chi 

Cm 



sin 

e 



sin^w 

e 



-\/2 sin 
e 



v^sin^ 



(7V2, + tan^wiVu) (C/,)*^ - V^YiN^, {Uj) 
2tan^wiV* (t/^)^„ + V2YiN^,^ {Uj)l^ 



Ra 



W 



(4.13a) 
(4.13b) 
(4.13c) 
(4.13d) 



Our results agree with those of refs. [11], [29]; the neutrahno contribution in ref.[12] seems to 
have some misprints. 

The analytic expressions of Sec. 2.2 can be used to rewrite both chargino contributions in 
terms of the loop functions /j and the basic SUSY parameters: 



SUSY 



me tan/3|/i|M2sin( 
"487r2sin2 9^ 



'=1 xf 



(4.14) 
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+ 2 [M|+ + 2tan/5M2|/i|cos0/, + 2M^cos2/3] J] ^^^^^"^ (^•^^) 
- \ [{Ml - - Y-^) + 2M^ cos 2/3(1 + 1^)] J] ^^^^ 

where Ac has been defined in eq.(2.11). Together with eqs.(4.12) these expressions exphcitly 
show that the chargino contributions to the leptonic dipole moments decouple like l/m|± for 
m|± ^ ml, and like 1/m? in the opposite limit m? S> ^|±- For completeness we have 
included terms oc Y^, even though eq.(3.14) shows that 1^ ^ 1; if these terms are neglected, 

(a^)gusY °^ ^ advertised earlier. 

Analogous statements also hold for the neutralino contributions, but because of the more 
complicated nature of neutralino mixing we were not able to find simple exact analytic expres- 
sions for these contributions. However, with the help of eqs.(2.17) and (2.18) one can derive an 
approximate expression for the neutralino loop contribution to dg. 



deY° e'm, f\Al + ^itanf3\ f2{mlJ\M,\') - Mm l/\M,\') 

' " •Sm(0i — 0e) 



e- J svsY SStt^ I cos^^lMil "^II ~ "^L 



tan / 7 sin (cpn + cpi) 
^ cos'^ew\Miii\ (|//|2 - |Mi|2) 



HM/2«/i^ir)- 



2 

/2«/|/.|^) 



tan/3sin0^ [\n 


'f2{ml/Ml) 


-MlMmlJlf, 




2sin^ /X ( A* 


2 - Ml) 



In the first line of eq.(4.16) wc have used an approximate treatment of selectron mixing, which is 
quite sufficient for the given purpose.* The last line in cq.(4.16), which involves the SU (2) gauge 
interactions, has a very similar structure as the chargino loop contribution (4.14); however, the 
overall factor in front of the neutralino contribution is four times smaller than that of the 
chargino contribution. Note also that eq.(4.16) does not contain contributions oc sin0i, which 
in our convention measures the relative phase between Mi and M2; only the phase of Mi relative 
to either the phase 0e in selectron mixing or to the phase of /i is relevant. 



*The electric dipole moment is chirality violating and hence proportional to the Yukawa coupling. Therefore 
slepton mixing, which is proportional to the Yukawa coupling, cannot be neglected. 
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4.3 Numerical analysis 

As well known [30, 31, 32, 11, 10], the experimental bound (4.1) on de provides stringent 
constraints on MSSM parameter space. For example, the chargino contribution (4.14) to 
can be estimated to be 



where msuSY stands for the relevant sparticle (sneutrino or chargino, whichever is heavier) mass 
scale. The chargino contribution by itself can therefore satisfy the experimental constraint (4.2) 
only for very small phase and/or very large sparticle masses. For sparticle masses not much 
above 100 GeV, one would need phases of order 10""^ (10~^) or less in the chargino (neutralino 
or slepton) mass matrices; if tan (3^1, this constraint would become even stronger. Such small 
phases are unlikely to lead to measurable effects in high-energy collider experiments [16, 17]. 
Alternatively one can postulate that sparticle masses are very large [32]. Since gaugino masses 
are coupled to parameters in the Higgs sector via one-loop renormalization group equations, 
whereas a similar coupling between first generation sfermion masses and the Higgs sector only 
exists at two-loop level [33] , naturalness arguments favor models with large slepton masses and 
relatively modest gaugino masses. The estimate (4.17) indicates that first generation slepton 
masses well above 1 TeV would be required if the relevant phases are 0{1). In that case these 
sleptons would be beyond the reach of the next linear e'^e~ collider, which will have center-of- 
mass energy TeV. Moreover, since FCNC constraints would then also indicate very large 

masses for second generation sleptons (recall that we assume them to be exactly degenerate 
with the first generation), a possible excess in could not be accommodated within the MSSM. 

We therefore focus on the third possibility for satisfying the constraint (4.2), where the 
different contributions to dg largely cancel [10, 11]; that is, the neutralino contribution must 
cancel the chargino contribution. Here we quantitatively analyze this possibility for three 
scenarios; later we will analyze high-energy observables that are sensitive to phases within the 
same scenarios. 

In all cases we assume that the ratio of M2 and \Mi\ is similar to that in models with 
gaugino mass unification at the GUT scale, which predicts [33] \Mi\ ~ O.5M2 at the weak scale. 
Similarly, we take values for the soft breaking masses of SU{2) singlet and doublet sleptons 
that are consistent with the assumption of universal scalar masses at the GUT scale. Recall 
that we assume degenerate first and second generation soft breaking parameters in the slepton 
sector: 




(4.17) 



men = ^Hr = ^la^ 

Ag — = A. 

The assumption of universal scalar masses at the GUT scale imphes [33] that 



(4.18a) 
(4.18b) 
(4.18c) 



mf ~mf +0.46M| 



(4.19) 
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at the weak scale. Finally, we are interested in scenarios where at least Ir, II, Xi X2 can 
be pair-produced at a "first stage" linear collider operating at ^/s — 500 GeV. 

This leads us to consider three different scenarios, which we call Bl, B2 and B3. Scenario 
Bl has lul = M2, i.e. is characterized by strong mixing between SU{2) gauginos and Higgsinos; 
this will occur in both the chargino and neutralino sector. Contrariwise, B2 has ^ M|, 
i.e. all Higgsino-gaugino mixing is suppressed. In these two cases we take a relatively large 
value of 1^41, which enhances slepton L — R mixing for small tan/?; we will see shortly that 
this increases the possibility of cancellations between the chargino and neutralino contributions 
to de- On the other hand, — mixing, while important for the calculation of dg, in all 
cases remains negligible as far as selectron production at high energies is concerned. Case B3, 
which is almost^ identical to the much-studied Snowmass "benchmark point SPSIA" [34], has 
intermediate gaugino-Higgsino mixing, as well as slightly reduced slepton masses. In all three 
cases we take four different values of tan (3. Moreover, we allow the three relevant phases 0i, 0^ 
and to float freely, i.e. we pick random values for these phases. (Recall that we work in 
a phase convention where M2 is real and positive.) These three scenarios are summarized in 
Table 1. Of course, we respect all relevant limits from direct searches for superparticles at 
coUiders, in particular at LEP2 [22]. 





\Mi\ 


M2 






1^1 




tan j3 


01, 0/., 0A 


Bl 


100 


200 


235 


180 


500 


200 


3, 6,9, 12 


e 


—71, 71 




B2 


100 


200 


235 


180 


500 


500 


3, 6, 9, 12 


e 


— TT, TT 




B3 


102.2 


191.8 


198.7 


138.2 


255.5 


343.2 


5, 10, 15, 20 


e 


— TT, TT 





Table 1: The three scenarios studied in this paper. All dimensionful parameters are in GeV. 

For simplicity and limited space of representation we only show results for two choices of 
tan/3 = 3 or 12 in scenarios Bl and B2, and for tan/3 = 10 or 20 in case B3. Results for the 
other cases are qualitatively similar and can be obtained by extrapolation from these extreme 
cases. Note that the numerical results shown below are projections of a three-dimensional 
parameter space onto two-dimensional planes. Hence it should be kept in mind that each 
correlation in the (p^ — (l)y plane has been obtained by scanning over the entire allowed range 
for (f)z. By far the strongest restriction on parameter space comes from de'. at least 99.4% of 
all randomly chosen points in a given run violate the constraint (4.2); the success rate at large 
tan /3 is even smaller.* In the following we will quote upper bounds on |0^| that result from the 
constraint (4.2). A similar band around 0^ = tt exists for small tan/3 and large |//|. 

Fig. 2 shows allowed combinations of the phases 0^ and 0i. We observe very strong con- 
straints on 0^ in scenario Bl, which become stronger as tan/3 increases. Scenario B2 allows 

'^The agreement becomes exact for the "benchmark value" tan/3 = 10 and vanishing phases. 

*This indicates that rather severe fine-tuning is required [17] to obtain the necessary cancellations if all 
phases are indeed independent quantities. To put it differently, one faces the challenge to construct models that 
"naturally" explain the required correlations between these phases. We will not attempt to do this here. 
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a: Bl, tan = 3 
1 



c: B3, tan/? = 10 

r 




-1 -0.5 0.5 1 

d: Bl, tan/5 = 12 
1 



1 -0.5 0.5 1 

e: B2, tan/3 = 12 

r 



-1 -0.5 0.5 1 
f: B3, tan/3 = 20 
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Figure 2: Combinations of 0^ and 0i that are allowed for at least one (pA £ [— ti", ti"]- 



much larger values of |0^|, which moreover do not decrease much with increasing tan/3, while 
scenario B3 is intermediate between these two. This behavior can be understood from eqs.(4.14) 
and (4.16). We saw that the contributions involving SU{2) gauge interactions have very sim- 
ilar structure in both cases, but the chargino loop contribution is bigger by a factor ~ 4 than 
this part of the neutralino contribution. The potentially most important cancellation there- 
fore occurs between the chargino contribution [more exactly: the total contribution involving 
SU{2) interactions, which is however always dominated by the chargino contribution] and the 
neutralino contributions involving U{1)y interactions. 

Scenario Bl has = M2, i.e. very strong mixing between Higgsinos and SU{2) gauginos. 
Eq.(4.16) no longer gives an accurate estimate of the neutralino contribution in this limit, 
but we expect it to remain qualitatively correct; note that it gives a finite answer (involving 
the derivative of the function in this case. In particular, the contributions involving the 
SU{2) gauge coupling would be much bigger than those involving U{1)y interactions if the 
relevant phases had similar magnitude; in other words, a significant cancellation can only 
occur if \(f)^\ is well below |(/)i|. Furthermore, for this choice of parameters a strong internal 
cancellation occurs between the two contributions from U{1)y interactions that grow oc tan/3, 
i.e. between the first line and the following two lines on the r.h.s. of eq.(4.16). As a result 
we find < 7r/30 even for tan/3 = 3. Moreover, the dominant contribution from f/(l)y 
interactions in this scenario involves A, i.e. is independent of tan /3, whereas the contribution 
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from SU{2) interactions increases oc tan/3. The upper bound on |0^| therefore scales essentially 
like cot (5. The importance of 0^ in this scenario also explains why there is almost no correlation 
between the allowed values of 4>^l and 0i. Moreover, in this scenario values of 0^ near tt are 
excluded by the lower bound (4.6) on (a^)susY- 

Eq.(4.14) shows that increasing \^\ while keeping all other parameters the same decreases 
the chargino contribution to d^] according to eq.(4.16) it also decreases the neutralino contri- 
butions that involve SU (2) interactions, but it actually increases the neutralino contribution 
that is sensitive to cl — mixing, i.e. the first line in eq.(4.16).^ Much larger values of 
\(f)n\ therefore now become possible. For the given choice of parameters the coefficient of the 
neutralino contribution oc sin(0i — 0e) is still about 5 times smaller than the coefficient of 
sin0^ in the chargino contribution, leading to an upper limit of ~ 7r/4 on Since these 

two coefficients have the same sign, cancellations obtain only if (pi + (p^ and (p^ have opposite 
signs. Note that both of these contributions are (essentially) oc tan/3. The upper bound on |0^| 
is therefore now almost independent of tan/3. However, one needs increasingly more perfect 
cancellations as tan/3 increases; moreover, the relative importance of the phase (pA diminishes 
with increasing tan/3, since its contribution to cl — mixing is not enhanced in this limit. 
These two considerations explain why the width of the allowed band decreases essentially like 
cot /3 for large tan /3. 

The increase of when going from scenario Bl to B2 also reduces the super symmetric 
contribution to a^. For tan/3 = 3 we therefore now also find an allowed band with (p^ ~ tt; 
however, this band disappears at tan/3 ~ 10. Note that the phase (pi enters mostly in the 
combination-i- cos(0i + (p^). This means that 0i ~ will give positive (negative) contributions 
to if 0^ ~ (tt). In other words, for values of (pi near zero the U{1)y interactions contributes 
with equal sign to as the (usually leading) SU{2) interactions do, whereas 0i ~ tt leads to a 
partial cancellation between U{1)y and SU{2) contributions. 0^ ~ tt therefore remains allowed 
to slightly higher values of tan /3 if 0i ~ tt as well. 

If \fj,\ is increased by another factor of ~ y/b, chargino and neutralino loop contributions to 
de can be of the same size, in which case no upper limit can be given on either |0i + or |(/)^| 
separately [11], although a strong (anti-) correlation between these two phases still has to hold. 
If is increased even further, the neutralino contribution becomes dominant. In that case (p^^ 
could take any value (after scanning over the other phases), but significant absolute constraints 
on the combination (pi + (p^ would emerge that hold even after scanning over all and (p^. 
However, most models of supersymmetry breaking prefer [33] values of |^| that are not much 
larger than M2. We therefore do not discuss scenarios with |//| ^ M2 any further. 

Scenario B3 has a significantly smaller value of than scenario B2, although it is larger 
than in Bl. The absolute upper bound on |0^| is therefore reduced to ~ 7r/8. The allowed 
bands in figs. 2c,f are narrower than in 2b,e due to the larger values of tan/3 and shghtly 
smaller slepton masses; both effects tend to increase the SUSY contributions to d^, requiring 

^In principle one can therefore have large cancellations between chargino and neutralino contributions even 
for M2 — 1^1, if M2 — m,y ^ \Mi\,mea- However, if M2 and are as in scenario Bl, this would require values 
of me„ well below the direct search limit of ~ 100 GeV. 

*For tan/3 » \ A\, cos((/>i — ^^) ~ cos((/>i + (^^) as well. 
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correspondingly more perfect cancellations. Note also that for tan/3 = 20 values of 01 near 
zero give above the range (4.6), i.e. in this case neutralino and chargino contributions to 
must not add constructively. Parameter sets with 0^ near vr are only allowed for tan /? ^ 5. 
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Figure 3: Combinations of 0^ and 0a that are allowed for at least one value of 0i G [— tt, tt]. 

The allowed regions in the (0^, 0^)— plane are shown in Fig. 3. Most of our scenarios have 
|yutan/3| significantly above 1^41, in which case the value of 0a is not very important. Even 



if 



')A is important, as in scenario Bl, there is little correlation between 0a and 0^, since 0a 



only enters in the combination 0e — 
on |( 



and 



1 is scanned in Fig. 3. In all cases the bound 
is slightly weaker for 0a — than for 0a — ±7r, since in the former case A and add 
(mostly) constructively to the mixing of selectrons, thereby increasing the first line on the r.h.s. 
of eq.(4.16). 

We see from figs. 2 and 3 that in all cases the entire range of values of 0a and 0i is allowed 
by the de constraint for some combinations of the other phases. Fig. 4 shows that there is little 
correlation between the allowed ranges of these two phases. Indeed, the df. constraint allows all 
combinations of these two phases, for some value of 0^. On the other hand, in case B3 with 
tan/? = 20 the constraint (4.6) excludes |0i| ^ 7r/2, see Fig. 2f. 

In Sec. 7.3 we will study correlations between low- and high-energy observables. To that 
end it is instructive to see how the low-energy observables correlate with the SUSY phases 
in the experimentally allowed region of parameter space. We saw above that 0^ is tightly 
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Figure 4: Combinations of 0i and 0a that are allowed for at least one value of 0^ G [— vr, vr] 



constrained, whereas 0^ and 0i are not. Since (pA does not affect high-energy observables, the 
most interesting correlations are those between low-energy observables and 0i, after scanning 
over (pA and 0^. 

We find that there is no correlation between de and 0i (not shown), whereas in scenarios 
B2 and B3 shows a behavior oc acos0i + b with a finite scatter, see Fig. 5. This difference 
originates from the requirement of very strong cancellations in discussed above. In particular, 
the phases 0i and 0^ have to be correlated such that the leading terms oc sin(0i + 0^) and 
oc sin 0^ cancel each other, to an accuracy determined by the size of (subleading) terms oc sin 0^ 
as well as by the experimental error on d^. This completely removes the correlation between d^ 
and sin 01 that one might naively expect from eq.(4.16). The phase-dependent neutralino loop 
contributions to (a^)susY can be read off from eq.(4.16) by replacing me by 2m^ in the overall 
factor, and all sin by cos; in addition, there are phase-independent contributions of comparable 
size. Since for our examples 0^ is constrained to be small (or near vr), | cos0^| ~ 1 and one 
finds a cos— like dependence of on 0i, as already stated. The crucial observation is that the 
01-dependent and 0^-dependent terms do not cancel in this case, so the 'naive' dependence of 

on 01 survives. We note in passing that (a^)susY = can usually not be achieved for a given 
choice of the absolute values of the SUSY parameters once we have required large cancellations 
in de, i.e. we cannot choose the phases such that there are large cancellations in both d^. and 
(o/i)susY- A better measurement of, and more accurate SM prediction for, therefore has 
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higher potential to further constrain the SUSY phases than improved measurements of de 
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Figure 5: (a^)susY vs. (pi after scanning over 0^ and (pA- 



5 High Energy Observables 

We are now ready to analyze the impact of SUSY phases on high energy observables. After 
defining the relevant kinematical quantities for the 2^2 processes under consideration, we 
briefly present the calculation of the corresponding unpolarized total cross sections. All these 
processes have already been discussed in the literature: results for e~e^ and e~e~ production 
results can be found in [35, 36, 37, 38] and [39], whereas results for XiXj ^^'^ Xi'x? production 
are given in [40, 36, 19] and [41]. We nevertheless list our results here in order to provide a 
self-contained presentation and to illustrate consistency with previous works. 



^ Of course, experimentally establishing a nonvanishing value of de would be of the greatest importance, since 
it would require physics beyond the SM. However, while it would require some SUSY phase to be non-zero, it 
would not further reduce the allowed ranges of any one of these phases after scanning over the other two phases. 
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5.1 Kinematics 
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Figure 6: Kinematical situation 



The kinematical situation is illustrated by Fig. 6. The momenta and helicities of the in- 
coming (first) electron and positron (second electron) are denoted by Pi and ui, and P2 and 
(72, respectively. The momenta of the produced superparticles, generically labeled by b and c, 
are denoted by k'^ and • case of fermions being produced their helicities are denoted by 
Ai and A2. 

Working in the center of mass (CMS) frame, we define the 2— axis of the coordinate system 
such that pi points in +z direction. The event plane is then completed by the momentum ki of 
particle b and defines the (x, z) plane of the coordinate system. The scattering angle 9 is defined 
as the angle between pi and ki. The nominal range for 6, which we use when going from the 
differential to the total cross section, extends from to vr. However, if the final state consists 
of two identical particles, physically 6 has to be < 7r/2; we therefore have to multiply the cross 
section for the production of identical particles with a factor of 1/2. Notice that our convention 
implies vanishing azimuthal angle 0. This definition of the {x, z) plane is convenient since we 
are only interested in total cross sections for unpolarized beams.* Of course, the phase space 
integration, which should be performed in a lab-fixed coordinate system, still gives a factor of 
27r from the integration over the azimuthal angle. Explicit expressions for the momenta and 
/cf can be found in Appendix A. 



5.2 Cross section for e e+ ^ e- et 

Fig. 7 shows the s- and t-channel diagram contributing to selectron pair production. By intro- 
ducing a dimensionless Z boson propagator 



*A nontrivial dependence on the azimuthal angle would arise only if we considered transversely polarized 
beams, and/or were interested in the kinematical distribution of the decay products of the produced superpar- 
ticles b and c. 
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Figure 7: Diagrams for e^e —>■ e 



7± 



and bilinear charges Z,^ 

(sin^ 6'w - 1)^ sin^ Ow - ^ ^ ^ 

ZJ, = l + 5ii^!^Z),. ZJ« = 1 + 4^C,. (5,2b) 

COS'^fc'w COS^fc'w 

^fi? = ^II = 0, (5.2c) 
the gauge contribution to the helicity amplitudes can be written as 

Ml''''''' = jv{p2,a2)Z^jP^u{p,,ai)ih - k.Y- (5.3) 
The neutralino contribution is 

= -V{P2, ^2)Kil,P^{^, - ^1 + m^o)DtKi,,P^u{pu a,). (5.4) 
The coefficients Kl^f. are given by 

^Lk = -7^ -— (cos^wA^2fe + sin^wiVifc), (5.5a) 

V 2 cos fc'w sm fc'w 

K^, = AT-,, (5.5b) 
V2 cos t^w 

K^k = Kk, = 0, (5.5c) 
and the neutralino propagators are 

(t, u) - mio 

where t = {pi — kiY and u = {pi — k2f' . By introducing a shorthand notation for the helicity 
amplitudes 

(a,a2).,=-Mr''' + -^r""' (5-7) 
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and using the explicit expressions for hclicity amplitudes given in Appendix B and the definition 
of the neutralino functions in Eqs.(C.2) we find six non- vanishing helicity amplitudes {9 is the 
angle between the momenta of the incident e~ and the produced e~): 



-2e'Ml^{s,ty, 



{++)rl 
{ )lr 

{+-)rr = -e'Xln sin 9 {Nrr{s, t) + Z+^) ; 
{-+)rr 



2e^MRL{s,t)- 



-e A|^sin^Z+^; 

1 

-e^Al^sin^Z^^; 



-e^\lLs\n9{NLL{s,t) + Zlj). 



(5.8a) 
(5.8b) 

(5.8c) 

(5.8d) 

(5.8e) 

(5.8f) 



Here, the kinematical factors = A|.g, are as in eqs.(A.3). As usual, the unpolarized cross 
section can be obtained by averaging over initial helicities. After integrating over the azimuthal 
angle, we obtain 







dcos9 


1287rs 


dCJRR _ 


1. 

^RR 


dcos9 


1287rs 


dcTLR 


1 

^LR 


dcos9 


1287rs 


daRL _ 


1 

^RL 


dcos9 


128ns 



H — )rr 
\{—)lr\' 
\{++)rl\' 



' + \{-+)RRn; 



(5.9a) 
(5.9b) 
(5.9c) 
(5.9d) 



Finally, for these and all following reactions the total, unpolarized cross sections may be ob- 
tained by performing the remaining integration over the scattering angle: 



dcos9 . 
I \ d cos 9 



da,j 



(5.10) 



5.3 Cross section for e e 



The t- and ti-channel diagrams contributing to pair production are shown in Fig. 8. The 
corresponding invariant amplitude can be written as 



= -^afe^Jfe^(P2,a2)x 



Pa + 5 



t—m 



IT- 



Pa N(pi,ai). 



(5.11) 
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+ 



Figure 8: Diagrams for e e e- 



Using tlie results of Appendix B we evaluate the helicity amplitudes and find 



(^1^2)., = --sin^^A55.,,_., {K.kKi^.kD^ - K^^^^kK^D^) 



(5.12) 



Rewriting these results in terms of neutralino functions as defined in Eqs.(C.2) we find four 
non-vanishing helicity amplitudes [9 is the angle between the momenta of an incident e~ and 
a produced e~; it does not matter which initial and final state particles are chosen, since the 
cross section is invariant under 6 —>■ rr — 6): 



{++)Rn = 2e' [M^^(s, t) + M^^(s, u)] ; 
( )ll = -2e^ [Mll{s, t) + Mll{s, u)] ; 

{-+)lr = e^xijiSmeNLR{s,t)- 

1 

(H — )rl = -e^XljiSm9NLR{s,u). 



(5.13a) 
(5.13b) 

(5.13c) 

(5.13d) 



After calculating the polarization averaged squared matrix elements and including the phase 
space factor, the differential cross sections are: 



dcTLL 




dcosO 


2567rs 


dcTRR 


_ ^RR 


dcosO 


256ns 


dcTLR 


1 

_ ^LR 


dcosO 


1287rs 



K 

{++)rr?] 

{\{-+)lr?+\{+-)rl?). 



(5.14a) 
(5.14b) 
(5.14c) 



Note that aiR and grl are not physically distinguishable in this case, unlike for e+e annihi- 
lation. 
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Figure 9: Diagrams for e ^ Xi X 



5.4 Cross section for e ^ Xi Xj 

Fig. 9 shows the s- and t-channel contributions to XiXj pair production. After a Fierz- 
rearrangement of the V contribution, the invariant amplitude can be written as 

-MS<x,;A,a, = ^v{p2,a2)i,P''u{p,,a^)Q%u,{hAi)YP^v,{hA2). (5.15) 
where we introduced the bihnear charges (with x = sin^ 9y^) 



gi^,^l + ^|^(.-^-lcos2,.), (5.16a) 



= 1 + , (5.16b) 

Qlr = 1 + ZtS,') - i - i -^20.) + ^(1 + cos 20,), (5.16c) 



ai , , / 3 1 



= 1 + YZr^ ( X - - - - cos20i ) ; (5.16d) 



i^z(2.T - 1) / 3 1 , , , 

,22 , , ^ 3 1 



= 1 + 1,"^ " i + i cos20,J , (5.17b) 

Dz(2x-1)/ 3 1 \ sD? , , , , 

Ql\ = 1 + -^^^^ - 4 + 4 + - (5.17c) 

,22 , Dz r 3 1 



^^^^ = 1 + ( X - - + - cos20^ ) ; (5.17d) 



Qll = {QlW = ^/[f ^^ in20,e-/^^ (5.18a) 
8x(l — xj 
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Qrr — {Qrr) 

Q^LR = {Q¥r) 

Qrl = (Qrl) 



sin 2(f) fie 



«(7l-/3jS-72) 



4(1 -x) 
Dzi2x - 1) sD 



8x(l — x) 
Dz 



Ax 



4x(l — a;^ 



sm zcpie 



(5.18b) 
(5.18c) 
(5.18d) 



Here the sneutrino propagator is defined analogously to the neutralino propagators (5.6). 
Using the results of Appendix B we find for a generic helicity amplitude [9 is the angle between 
the momenta of the incident e~ and the produced x~): 



(5.19) 



+ h^-\2\l (1 + /5Air7^Ai)(l + f3Xir]_f3Xi) {cosO + Aidi) 



where the kinematical quantities rj± are defined in eq.(B.3). The unpolarized cross sections can 
be computed from eq.(5.19) by averaging over initial helicities and summing over the final ones: 



daij 



Txa' 



d cos 8^ -^1 { [(1 - ) + cos^ 6] + 8/i./i,Q^ + 2a| cos OQ^ } , (5.20) 
where 47ra = e^, /ij = m^J ^/s and Ajj = nf — ji^-. The new quartic charges are given by: 

Qi' = \Qt+? + \Q%-? + + \Q'--?; (5.21a) 

Q^^' = {QX+Qf_ + Q'1.Q%) ; (5.21b) 

Q'i = \QW - \Qt\' - IQ^P + \QW- (5.21c) 
5.5 Cross section for e e+ X^X^ 




+ 



Figure 10: Diagrams for e e"^ x^X^ 



In Fig. 10 the s— and t— channel contributions to XiXj production are shown; the addi- 



tional, destructively interfering -u— channel diagram is indicated by the exchanged indices in 
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parentheses. Applying a Fierz rearrangement on both the t— and m— channel diagram, and 
re-ordering the channel amplitude, we obtain the invariant amplitude 



A^S.;A,A. = ^^(P2, '72h,PMPi, <^2)Q%u{k,, Xih'P'^vih, X,). (5.22) 
Here, the bilinear charges are given by {x — sin^ ^w) 

= -;wT^(2^ - 1)% + ^^t9up (5-23c) 



2x(l - x) 



= ^^i;- + sD^gn,,, (5.23d) 



with 



to,- = ^ (A^2* + tan ^wiVii) (iV2i + tan ^w^Vi,) , (5.24a) 

The selectron propagators are defined as 

Dti!" = / 2 • (5-25) 

Since this amplitude has the same structure as the amplitude for XiXj production, eq. (5.15), 
we can directly translate the result (5.19) from this calculation; we just have to replace the 
bihnear charges. We can also use the result (5.20) for the unpolarized, differential cross section, 
but we have to include a statistical factor: 



d cos 9 8s 



a| I [(1 - Aj) + Ai,- cos=^ e] + ^[ii[ijQi + 2a| cos QQ'i } . (5.26) 



Of course, now the bilinear charges of eqs.(5.23) have to be used when evaluating the quartic 
charges defined in eqs.(5.21). 



5.6 Pat for 2 fermion production 

We will see in Sec. 7 that some of the cross sections calculated in the previous subsections depend 

quite sensitively on the CP-violating phases (/>! and/or 0^. Nevertheless, if measurements 
of these cross sections establish a deviation from the CP-conserving MSSM, one will need 
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to measure some CP-violating asymmetries in order to convince oneself that the observed 
deviation is indeed due to non-vanishing phases, rather than due to some extension of the 
MSSM. We will see that this is possible only for the production of fermionic final states. 
Consider 

e"(pi, Si)e"^(p2, S2) Xi{ki, h)Xj{k2, h) ■ (5.27) 

— * — * 

The momenta pi_2 and ki^2 have been defined in Fig. 6, and Si^2 a-nd Si^2 are the spin vectors 
in the initial and final state, respectively. A CP transformation on reaction (5.27) gives the 
CP-conjugate process 

CP : e+(-pi, si)e~(-p2, ^2) Xi(-ki, li)Xj(-^2, ^2) ■ (5.28) 

In the center-of-mass system, pi = —p2 and ki = —k2- The initial state will therefore be self- 
conjugate if si = S2, in particular for unpolarized beams. Comparing reactions (5.27) and (5.28) 
one can introduce two CP-odd asymmetries even after summing over the spins in the final state. 
One can define a rate asymmetry for chargino production, essentially cr{xixt) ~ "^(x^Xi")) as 
well as an angular asymmetry for the production of two different neutralinos, proportional 
to da{x^x^,6) — da{x^x^,TT — 9). However, far from the Z pole, both these asymmetries 

vanish identically at the tree level. The reason is that they are odd under a combined CPT 
transformation, where the "naive time reversal" T reverses the direction of all 3-momenta, but 
does not exchange initial and final state. Quantities that are odd under CPT can be non-zero 
only in the presence of absorptive phases, which can come from nearly resonant s— channel 
propagators, or from loop corrections if the kinematics allows the particles in the loop to be 
on-shell. 

A CP-odd quantity can therefore only be non-zero in the absence of absorptive phases, if it 
is also T— odd. This is true for triple products of momentum and spin vectors. In general, the 
spin of the final state fermions in (5.27) can be decomposed in three components: P^'*-' is the 
component of Sj in direction of fcj, averaged over many events (with fixed 9)\ P^'*-' is orthogonal 
to fcj, but lies in the event plane; and Pjy*"' is orthogonal to ki anc? orthogonal to the event plane. 
The first two of these quantities arc T— even; however, since 

= (si-(plx^,)), i = l,2, (5.29) 

P^*-' is indeed T— odd; here (. . .) denotes averaging over many events with fixed scattering angle 
6.* We will comment on the measur ability of this quantity when we present numerical results. 

The normal components of the polarizations of Xi and Xj can be computed using results of 
ref.[42]: 

_ -2^ni {E.,., [(^1^2; + + ) ((71(72; - + Y + {(Jl(J2] +') {(Jl(J2] —)"] ] _ . 

|/„ „ . , , \ 2 I „ . , \I2 , l/„ ^ . , \I2 I 1/ . r|2' lO-iUaj 



E.,.J(^i^2;++)P+ 1(^1^2; + \{a,a2; + \{a,a2; 

* Strictly speaking, is CP-odd only for self-conjugate final states (any two neutralinos, or Xixl). How- 
ever, since at tree-level and away from s— channel resonances T and T transformations are essentially the same, 
a non-vanishing in XiXt production can also be considered evidence for CP-violation. 
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E.,.. [1(^1^2; + 1(^1^2;+-) |2 + |(aia2;-+)P + |(aia2;—)P]- ^ " ^ 
After introducing a fourth quartic charge, 

Q^' = {QX+Qf_ + Q'LQ%) , (5.31) 

and using eq.(5.19), eqs.(5.30) simphfy to 

P'^'' = ''^^ ; (5.32a) 

[1 - A?. + X,j cos2 9] QY + SmujQ'i + 2X1 

Pi^^ = , = -P^^^. (5.32b) 

[1 - A2. + \j cos2 ^] + ^li^iijQ'i + 2A[;. COS ^^Qg^ 

We see that P^"' vanishes both at threshold (where Ajj — 0) and far above threshold (where 
~^ 0). Eqs.(5.16) and (5.17) show that all bilinear charges describing XiXt production are 
real. Likewise, eqs.(5.2) and (5.24) show that the couplings appearing in the expressions (5.23) 
of the bilinear charges for neutralino pair production are real for final states consisting of two 
identical neutrahnos. We thus see that Q^l and hence P^-'^''-' can only be non-vanishing for 
off-diagonal production modes {i 7^ j). Moreover, the second identity in eq. (5.32b) shows that 
there is only one independent P/v for each distinct XiXj production channel, for a total of 7 
independent CP-odd observables. 



5.7 Approximate results 

The results presented in the previous sections allow the exact (tree-level) calculation of the 
phase dependences of the selectron and neutralino production cross sections, and of P^r. How- 
ever, it is useful to get a qualitative understanding of where one can expect strong sensitivity 
to the fundamental phases in the supersymmctric Lagrangian. To this end we here discuss 
the behavior of the relevant cross sections and polarization components using the approximate 
diagonalization of the neutralino mass matrix described by eqs.(2.17)-(2.19). 

We begin with the cross sections for selectron pair production. All modes receive (9(M^) 
contributions from the exchange of the bino-like neutralino; in case of cl pair production the 
exchange of the wino-like neutralino also contributes at order M^. The e^e^ mode is the only 
one which has a phase sensitivity to order M^, where the cross section is sensitive to the relative 
phase between Mi and M2. The other cross sections show phase sensitivity only at order M|. 
This is due to both the exchange of heavier, Higgsino-like neutralinos which develop gaugino 
components at 0{Mz)., and due to the C(M|) corrections to the gaugino components of the 
gaugino-like neutralinos; the importance of these latter contributions explains why we included 
the (9(M|) quantities 812 and ^211 arid included the normalization factors A^i,2) in eqs.(2.17) 
and (2.18). The 0{Mz) shifts 5mi^2 of the masses of the gaugino-like neutralinos also affect the 
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selectron production cross sections, either directly (if the physical masses are allowed to vary 
with the phases), or indirectly (if physical masses are kept fixed, in which case the absolute 
values of the input parameters have to be varied along with the phases); since 0(M|) shifts 
of |Afi|, M2 and change the eigenstates only at (9(M|), we will ignore such indirect 
effects in the following. Since does not have SU{2) interactions, o"(e^e^) are at 0{M'^) 
only sensitive to the phase combination 0i + 0^, whereas the other modes are also sensitive 
to* 01 and 0^. One should also bear in mind that the phase sensitivity of the diagonal e'^e~ 
production channels is further diluted by the presence of large s— channel (7 and Z exchange) 
contributions, which do not depend on any phase. 





COS(0^ + 01 ) 


COS0^ 


cos 01 




sin 2/9^^1*^^! 










M2IMP 


^L^R 






M|M2 
IMiIImP 


^R^R 


sin2/?-flfl 


none 


none 


^L^L 






|Mi| 
M2 


^L^R 






M||Mi| 








^R^R 




none 


none 



Table 2: Phase dependence of the cross sections for selectron pair production in e+e" as well as 
e~e~ annihilation for fixed physical neutralino masses. Each entry gives the dependence of the 
coefficient of the indicated (combination of) phase(s) on the supersymmetric parameters relative 
to the leading (phase-independent) contribution to this cross section, under the assumption 
|Afip < M| ^ l/ip. "None" means that the corresponding term does not exist to (!?(Af|). The 
cross section for e^e^ production also has terms cx cos(20i + 0^) and oc cos(0^ — 0i), but with 
small coefficients oc sin 2/3Af|/(|/i|^). 

The phase dependence of the selectron production cross sections for fixed physical neutralino 
masses is summarized in Table 2. Here we show the coefficients of the various phase-dependent 
terms that can appear, relative to the leading (phase-independent) contribution to this cross 
section. We have omitted numerical factors, including factors involving the weak mixing angle. 
Nevertheless we can draw some conclusions from this table. First, we notice that the dependence 
on the phase 0^ shown in the second and third columns vanishes like 1/ tan/? for tan (3 ^ 1. 
The reason is that the dependence on this phase in the neutralino mass matrix could be rotated^ 

*Recall that 02 = in our convention. 

^This rotation does not introduce any phase in those parts of ffx vertices that come from gauge interactions, 
but does introduce a phase in the Yukawa contribution to these vortices. Recall that these Yukawa contributions 
can be ignored when calculating cross sections, but have to be kept when computing leptonic dipole moments. 
This explains why the 0^ dependence of dg and ttfj, is not suppressed at large tan/3. 
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into the off-diagonal gaugino-Higgsino mixing entries oc cos /?. However, tlie dependence on 
the relative phase between the two soft gaugino masses does not vary with tan (5. 

Second, with the exception of the e£e£ mode, all phase dependence vanishes as |//| — > 
oo, but the dependence varies for different modes. In particular, the phase dependence 
of the diagonal mode e^e^ vanishes oc for large whereas all other cross sections 

receive phase-dependent contributions that only fall like however, for tan/3 — > oo the 

I dependence of the total phase sensitivity becomes stronger, as can be seen in the last 
column. In most cases the leading phase dependence comes from the exchange of the lighter, 
gaugino-like, neutralinos. The exception is the e^e^ mode, where the the exchange of the 
heavier, Higgsino-like states contributes at the same order. 

Clearly the LL mode will have the strongest phase dependence of all e~e~ channels [38], 
and indeed of all selectron production channels, since it already occurs at 0{Mz), as noted 
earlier. The phase-dependent terms in e£e^ and e£ej production are of similar size. For our 
choice of parameters the second mode is preferable, since it is accessible at lower energies, and 
since the cross section near threshold scales like V^, rather than like A^^^. Finally, the relative 
importance of phase-sensitive and phase-insensitive terms in most selectron production cross 
sections does not depend strongly on the beam energy. We therefore expect the best statistical 
accuracy for the determination of the relevant phases when the beam energy is chosen such 
that the cross section being investigated is maximal. 

So far we have kept the physical neutralino masses fixed, which means that \Mi\, M2 and 
1/^1 have to be varied along with the phases; we saw above that this affects the cross sections 
only at C(M|) relative to the leading term. If instead these input parameters are held fixed, 
the physical neutralino masses will vary at 0(M|). Of particular interest are the masses of the 
gaugino-like states, whose exchange gives much bigger contributions to the matrix elements 
than that of the Higgsino-likc neutralinos. The relevant mass shifts are given in eqs.(2.19). We 
see that these effects also vanish oc 1/ tan/3 for large tan/3. However, they only scale like 
for large They will therefore dominate the total phase dependence of the e~^e]^ production 
cross section. For the other modes, the dependence on cos(0i + 0^) and on cos0^ that comes 
from the variation of the masses of the gaugino-like neutralinos is qualitatively the same as 
shown in Table 2, if we ignore factors \Mi\/M2. A more detailed analysis is therefore required 
to decide which source of phase dependence dominates. However, in case of ej^e]^ production 
the total phase dependence is still dominated by the 0{Mz) term from bino-wino interference. 

We now turn to the cross sections for neutralino pair production in e'^e~ annihilation, aij = 
a{e^e^ Xi'x^)- We first note that of the 10 distinct cross sections, only four receive 0{Mz) 
contributions: the cross sections an, 0"i2 and (T22 describing the production of two gaugino- 
like neutralinos receive large contributions from selectron exchange in the t— or m— channel, 
while (734 receives large contributions from Z exchange in the s— channel. The cross sections 
(733 and (T44 describing the production of two equal Higgsino-like states receive non-vanishing 
contributions only at C(Af|), whereas the cross sections for the production of one Higgsino-like 
and one gaugino-like state start at (9(M|). 

Only (J12 has sensitivity to some phase (in this case, 0i) at order . All other cross 
sections are sensitive to phases only at order M| or even M|. The strong phase sensitivity of 
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ai2 can be traced to the (^2— term in eq.(5.26). It comes from the fact [31] that the production 
of two Major ana fermions is P— wave suppressed near threshold if they have the same relative 
CP-phase, whereas any difference in this phase leads to an 5"— wave contribution to the cross 
section. This effect can be probed with optimal statistical significance rather close to threshold, 
in this case for ^/s not too much above \Mi \ + M2. 

The 0{M'^) phase-dependent terms in the neutralino production cross sections should be 
most easily observable in the mixed "gaugino-Higgsino" final states, since here the cross sections 
also only start at 0(M|), as remarked above. Note that the two Higgsino-hke neutrahnos are 
closely mass-degenerate in the limit 3> Mf . This makes it very difficult to experimentally 
distinguish between the production of X3 and x% In the following discussion we therefore always 
sum over these two Higgsino-like states. Once this has been done, we again find that all terms 
involving cp^ come with a factor sin 2/3, and are thus suppressed at large tan/3. These cross 
sections also contain terms oc cos(0i — (p/^) and cos(20i + 0^), which result from the rephasing- 
invariant combinations of phases ±(0i — (j)2) — {<t>2,i + <t>fi) in our convention 02 = 0. Altogether 
we find the following phase-dependent terms in these two cross sections: 

^ih = ^13 + (^14 : cos(0^ + 0i) I , j ; COS0, 



sin 2/]p/i 






li\s 



cos(0i - 0^) ; cos 01 I ,„ , ; (5.33a) 



\IjL\S \ S 



\Mi\ 


M2 






|2 



^ , /sin2/3|Mi| sin2/3|Mi|M|\ , /sin2/3M2 sin2/3|/i|M2\ 

^2h = ^23 + 0-24 : COS(0^ + 0i) I r- , p ; COS 0^ ( , ; 

sin2^M^ f \M,\M2 \M1\M2\ 
cos(20i + 0^) ; cos 01 I ,„ , . (5.33b) 



A common factor oc a^Mf /(|/xps), characterizing the size of the leading phase-independent 
contributions, has been factored out. Here we have listed the s— dependent contributions com- 
ing from the terms oc Q2 separately, where present. Note that they usually have a different 
dependence-'- on \fi\ than the terms that survive for s 00. We conclude from eqs.(5.33) that 
(7^^ might show a somewhat stronger overall phase dependence in the region of parameter space 
allowed by low-energy data, since it depends on the potentially large phase 0i + 0^ through 
terms with fewer powers of |//| in the denominator than a"2^ does, whereas the dependence on 
cos 01 is parametrically the same in both cases. This is fortunate, since the Ih mode is accessi- 
ble at lower energies. Finally, note that the phase dependence of the neutralino masses affects 
these cross sections only at 0{M^). 

We do not list the C(M|) phase dependent terms of the cross sections that receive C(M|) 
contributions, since these will clearly be much more difficult to measure. 

The situation concerning chargino pair production is rather similar. Here both diagonal 
modes start at 0{Mz), but receive phase-dependent contributions only at 0{Mz). In case 

tin some cases these threshold terms seem to grow with increasing |^|. However, (t^^ {'^2h) '^^ accessible only 
for y/s > |Mi| + ImI (•\/s > M2 + |/Lt|), i.e l/xp/s < 1 in the physical region. 
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of the ofF-diagonal mode^ both the cross section and the phase dependence starts at 0{M'^); 
indeed, the phase dependence is very similar to that in eq.(5.33b) with 0i — > 0, since the ^7(1)^ 
gaugino mass does not appear in the chargino mass matrix. 

Finally the results for perturbative neutralino mixing may be applied to the polarization 
vector components of the neutralinos produced in e~^e~ Xi'Xj calculated in Sec. 5.6. Here 
we only discuss the normal component as it is the only CP-odd quantity available if neutralino 
decays are not included explicitly. Recall that a non-vanishing P/v can only occur for final 
states consisting of two dijferent neutralmos. We find that the numerators in eqs.(5.32) receive 
0{M^) contributions only for the (12) mode; in case of the (l/i) and (2/i) modes the numerator 
starts at (9(M|), just like the corresponding total cross sections, and hence the denominators 
in (5.32). In all these cases P/v will therefore receive 0{M^) contributions. On the other hand, 
Pn for the hh [or (34)] mode vanishes to 0{Mz); this final state is therefore of little interest 
in the present context. Exphcitly, for (12) production we find to 0(M^):^ 



sm ( 



\^ - ^12 + ) J^LJ^ 




where z = cos 9. As expected for a CP-odd quantity the dominant dependence on (pi is through 
a sine function, while the denominator (basically the differential cross section discussed above) 
contains a CP-even dependence on 0i through a cosine. 

Due to their mass degeneracy we have to average P/v for the mixed gaugino-Higgsino modes 
over the production of both Higgsino~like neutralinos. Using the event numbers N as weights, 
we obtain: 

" " It;'^ , ^=1,2. (5.35) 

iVj3 + iVj4 

This amounts to replacing the quartic charges in eq. (5.32) by: 

Qt = Qf + Ql'- (5.36) 

The calculation of the relevant quartic charges Qf' to 0{M^) is now straightforward, if some- 
what tedious. We find the following terms, factoring out Q;^M|/|/ip: 

^-.r . n sm2p\Mi\M2 , ,sin2/5M| M2 , 

Q4 : sm(0i + 0^) sm2/i; sm0^ — ; sm(0i - 0^) — — ; sm0i— — ;(5.37a) 

^2h /, , ,sin2/5|Mi|M2 . , . . ^, ,sin2/5|MiP 
Q4 : sm(0i + 0^) — ; sm(p^sm2(3; sm(0^ + 20i) — ; 

1/^1 |A*I 

sin 01 ^^]^ . (5.37b) 
1^1 



Hn principle X1X2 production is now distinguishable from X1X2 production. However, the two cross sections 
differ only in the presence of an absorptive phase, i.e. after including loop corrections. 

^Eqs.(5.32) show that \Pj^ \ is larger than \P^ \ by a factor M2/IM1I ~ 2. However, we assume that Xi is 
the LSP, and hence stable (if R-parity is conserved), so that its spin cannot be measured. 
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The terms in eqs.(5.37) directly correspond to terms in P^-, up to an additional factor of 

\Mi\/\f^ (l/^l/v^) foi' P'}f^'^ {Pn^^)i since the dependence of the leading term in the denomi- 
nator in eq.(5.32) on SUSY parameters has already been factored out. As expected, the phase 
dependence is through sine functions here, and all terms that are sensitive to 0^ are suppressed 

at large tan j3. The first term in QY^ gives rise to a contribution to P^^^ that remains finite as 
— > oo, but vanishes oc 1/ tan/9 for large tan/5. On the other hand, in the {2h) mode we can 
measure the polarization of the lighter gaugino-like neutralino, giving rise to an extra factor 
\lJi\/\fs. We thus see that (as long as y/s > \fi\ + M2) the second term in eq.(5.37b) gives a 

contribution to P^'^'* that rises with increasing However, for the range of of interest to 
us, this term is suppressed by the stringent upper limit on | sin0^|, see Sec. 4.3; only in scenario 
B2 with small tan/3 can it reach comparable magnitude as the last term in eq.(5.37b). This 

last term leads to a contribution to P^'^'^ that approaches a constant for large ancZ remains 
finite for large tan/5. We therefore conclude that {2h) production should allow a somewhat 
more sensitive direct probe of CP violation than (Ih) production. Finally, the normal compo- 
nents of the polarization vectors in X1X2 production have similar structure as eq.(5.37b) with 
01 — > 0, but receive additional contributions from the Z coupling to the gaugino component of 
the heavy chargino state xf- 



6 Significances 

Oiu aim in this Section is to introduce objects quantifying the impact of CP-odd phases on 
total cross sections, which are CP-even quantities. To this end we compare the difference in 
counting rates between a CP-conserving point in parameter space (CPC: all phases (pi — 
or tt) and a CP-violating one (CPV: identical absolute values, but (pi ^ and low-energy 
compatible) to the statistical error at the CPC point. This determines the significance S with 
which a deviation from the cross section predicted for the CPC point can be measured. It can 
be written as 

^ _ AiVcpc-cpv _ A^cpc ~ Ncpy 
5Ncpc "//Vcpc 

Since there are two CP-conserving values (0, tt) for each phase, we have to deal with eight CPC 
points for each set of absolute values, and hence the same number of significances is available 
for each kinematical accessible cross section. The smallest of these evidently determines the 
statistical significance with which the presence of CP-violating phases can be inferred from 
this cross section for given values of the absolute values of all SUSY parameters. We therefore 
define as our final measure of the sensitivity of a given cross section to phases the significance 

|(jCPV _ ^CPC„ 

S (fifj) = min„ ( ^'^'r^^' ' | X VZ, (6.2) 
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where cx/./^ is the total cross section for e~e^ — > fifj, and n = 1, . . . , 8; we only include CPC 
points which are low-energy compatible.^ Finally, C is the integrated luminosity, which is 
expected to be different for the e+e~ and e~e~ options. 

In the procedure outlined so far, the CPC and CPV points have the same absolute values of 
Ml, M2 and fi. This means that these points will in general have different physical neutralino 
and chargino masses [43]. Recall that the phase dependence of the x masses is suppressed by 
Mz/{\fJ^\m^), see eq.(2.19). Nevertheless, changes of several percent are possible, in particular 
in the neutralino sector. This could lead to similar changes in the cross sections through kine- 
matical factors (in x production) or through neutralino propagator factors (in e production). 
Moreover, these masses are often more easily measurable than the cross sections which are the 
focus of this analysis. 

We therefore introduce a second set of significances S where CPC and CPV points have 
the same physical masses for XitXs ^iid xf'i m the limit of large x masses and for our choice 

> M2 > \Mi\, these three masses essentially fix |Afi|, \fi\ and M2, respectively. Note that 
we only have three (dimensionful) absolute values that can be adjusted in the neutralino and 
chargino mass matrices. We can therefore not guarantee that all chargino and neutralino masses 
are the same in the CPC and CPV points. However, after ensuring that these three x masses 
are the same in both points, the remaining variation of the other three x masses between the 
CPC and CPV points is quite small. For technical reasons we keep \Mi\, M2 and |/i| fixed (at 
the values listed in Table 1) for the CPC points, and adjust them at the CPV points. Since the 
eight CPV points have four different x mass spectra, a given set of phases now also produces 
several different CPC points. The new significance can thus be written as 



I CPV„ _ ^CPCnl \ 

S (/,/,) = min„ I r^^' >< (6-3) 

^fd, " / 

Our algorithm for calculating the significances can be summarized as follows: 

- Select a CPV point. For a set of the absolute values of the relevant SUSY parameters, 
as listed in Table 1 for our three scenarios Bl, B2 and B3, this amounts to randomly 
choosing values for the phases (pA, 4'n ^-^id 0i. Repeat this step until a point that is 
compatible with the low-energy constraints has been found. 

- For each process, find the low-energy allowed CPC point that minimizes S{fifj) as defined 
in eq.(6.2). Note that there are only eight CPC points for each scenario Bl, B2 and B3 
if tan P is kept fixed; however, this procedure in general selects different CPC points for 
different processes. This completes the calculation of the S. 



Define four new CPV points CPV^ by adjusting \Mi\, M2 and such that m^o, m^± 
and m^o are the same in points CPV„ and CPC„. 



^Since Uf^f. does not depend on tp^, there are only four different values of f^.j*^" for a given CPV point. 
However, occasionally both <j>A = and (pA = have to be checked to find a CPC point that is compatible with 
the bound on a^. Of course, the bound on dg is trivially satisfied by all CPC points. 
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- Calculate the S (fifj) as in eq.(6.3). 

Note that S and S only measure statistical significances. In addition there will be systematic 
uncertainties, both from experiment and theory. We have little to say about experimental 
systematic errors, except that we hope that they will be small. A theoretical error is introduced 
since our cross sections can only be predicted with finite precision. At tree-level these cross 
sections are determined uniquely by the parameters listed in Table 1 , plus a few SM parameters 
that are already now known with high precision. However, explicit calculations for xf pair 
production show that quantum corrections can easily amount to (9(10%) [44]. Some of these 
corrections can be calculated unambiguously once the parameters listed in Table 1 are specified, 
but the remaining corrections can still amount to several percent. In particular, the lepton- 
slepton-gaugino "gauge couplings" depend (logarithmically) on the squark mass scale [45]. 
The production of Higgsino-like charginos [44] and, presumably, neutralinos also depends on 
the parameters appearing in third generation sfermion masses. These corrections will only 
be calculable once the parameters of the (presumably quite heavy) squark sector have been 
determined. Until this has happened, out of two processes with roughly equal significances as 
defined above, the process with a smaller cross section should be preferred, since here a given 
significance corresponds to a larger relative variation of the cross section with the phases. 

7 Numerical Analysis 

We are now ready to present numerical results for our high-energy observables. We will first 
discuss the impact of the CP-phases on the (CP-even) cross sections, before turning to the (T- 
odd) normal components of x polarization vectors. Finally, in Sec. 7.3 we will study correlations 
between phase-sensitive quantities. 

7.1 Cross sections 

As discussed in Sec. 4.3, we chose our SUSY parameters such that selectron pair production 
as well as the production of two lighter neutralinos or charginos is possible already at the first 
stage of a future linear e+e~ collider (LC) operating at ^/s = 500 GeV, which is our default 
choice. However, in scenario B2 the Higgsino-like states are not accessible at this energy. In 
this scenario we therefore take y/s — 800 GeV when discussing reactions where at least one 
X3) X4 or xt state is produced; note that all current LC designs foresee an upgrade to at least 
that energy. A similar treatment is used in scenario B3, except for the X1X34 final state, which 
is already accessible at ^/s = 500 GeV in this case. 

In Table 3 we show the maximal allowed cross sections for the 19 different production 
channels discussed in Sec. 5, for our three scenarios Bl, B2 and B3, and the same choices of 
tan P employed in Sec. 4.3. Only combinations of phases that are allowed by the low-energy 
constraints on and have been included in the maximization. These cross sections have 
been calculated at tree-level, as described in Sees. 5.2-5.5. We have also ignored corrections due 
to initial-state radiation and beamstrahlung. These effects are often larger than the dependence 
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on CP-violating phases; they should therefore certainly be included in any future experimental 
analysis (along with radiative corrections, which will likely be known well before the first LC 
commences operations). However, they are largely independent of CP-phases, and should 
therefore not affect our conclusions. 
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Table 3: Maximal values of the total cross sections [in fb] for unpolarized beams, for the 
scenarios defined in Table 1. "-" means that the corresponding mode is not accessible. In sce- 
narios B2 and B3 we have summed over the production of the heavy Higgsino-like neutralinos, 
as described in the text. The beam energy is 500 GeV in most cases, but has been raised to 
800 GeV for the production of X3 4 states in scenarios B2 and B3, as indicated by the 

asterisk. Note that the charge-conjugate mode is included, if it is distinct from the listed one. 

We saw in Sees. 2.3 and 5.7 that the two heaviest, Higgsino-like neutralinos are close 
in mass if |/^| > M2 and ^ M|; the degeneracy between these states is only lifted at 
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0(M|/[|/ip — M|]) (as well as by radiative corrections, which however are sizable only in the 
presence of large A— terms in the stop sector [46]). Numerically, we find that the relative 
difference between m^o and m^o ranges from 24 to 35% in scenario Bl, but only from 0.2 to 
3.5% (0.1 to 7.5%) in B2 (B3). Since the production of nearly degenerate particles is difficult 
to distinguish experimentally, we simply sum over the production of Xs and X4 in scenarios B2 
and B3; in particular, we only give results for a single process of heavy Higgsino-like neutralino 
pair production in these cases. Recall that we used the same treatment in Sec. 5.7, eqs.(5.33) 
and (5.35)-(5.37). 

As well known [39, 41] many of our cross sections can be enhanced by factors of a few 
if both beams are polarized. Moreover, the discussion of Sec. 5.7 indicates that the greatest 
sensitivity to phases comes (through 0i) from the interference of SU{2) and U{1)y interactions; 
these contributions will be suppressed if one chooses beams, since is a singlet under 
SU{2). However, the sensitivity to other combinations of phases is enhanced for different 
choices of beam polarizations. We therefore only show results for unpolarized beams, with the 
understanding that in many cases the cross section (phase sensitivity) could be enhanced by 
up to a factor of 4 (2) if fully polarized beams were available. 

We see from Table 3 that the cross sections for selectron pair production are generically 
bigger at e~e~ coUiders than at e'^e~ colliders [37]. This difference is only partially compensated 
by the higher e~^e~ luminosity; we assume J Cdt — 500 (100) fb~^ for e^e~ (e~e~) collisions. 
We use these relatively conservative values since we do not include efficiency factors. These are 
expected to reduce the actually available event samples by factors of a few, the precise values 
depending on both the process under consideration and the sparticle spectrum. Moreover, at 
e'e" colliders the diagonal, chirality-conserving modes have higher cross section than the off- 
diagonal, chirality-violating mode; recall that the latter is P— wave suppressed near threshold, 
and vanishes for vanishing gaugino masses. At e+e~ colliders the diagonal selectron production 
modes are P— wave suppressed; this explains the rather small cross sections for e^e^ production. 
Finally, the selectron production cross sections are highest in scenario B3, since the selectron 
masses are somewhat smaller than in the other two cases; this effect is particularly significant 
for e^e^ production, which is a P^wave process quite close to threshold. The strong tan/3 
dependence of the maximal e^e^ production cross section in this scenario follows from the fact 
that the region near 0i = 0^ = is excluded by the constraint for tan/3 = 20, see Fig. 2f. 

The biggest cross sections at e+e" coUisions are those for e^e^, XiXi and x?x? production. 
However, the latter leads to an invisible, and hence undetectable, final state if is a stable 
LSP; we will therefore not analyze it any further. The cross sections for producing two heavy 
charginos or neutralinos are suppressed both by phase space and by their Higgsino-like nature. 
However, the production of one light and one heavy x state is possible in all three cases. Since, 
as discussed in Sec. 5.7, these cross sections are non-vanishing only in the presence of gaugino- 
Higgsino mixing, they fall with increasing |//|. However, even in scenario B2 one will have 
several thousand events containing these Higgsino-like states. For the other channels, typically 
several tens of thousands of events will be available, meaning that the cross sections could be 
measured with statistical uncertainty of 1% or less. 

The maximal possible values of the significances S and <S of eqs.(6.2) and (6.3) that can be 
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found in our three scenarios are summarized in Table 4. The e^e^ mode shows the strongest 
phase dependence of all selectron production channels, i.e. the highest significance, largely 
independent of |//| and tan/3; the tan/3 dependence of S in scenario B3 is due to the fact that 
the point </>! = (/)^ = is excluded by the constraint at tan/3 = 20, but still allowed at 
tan/3 = 10, as shown in Fig. 2c. The mixed e^e^ mode is the for our purposes most promising 
selectron production mode at e"'"e~ colliders. It would allow to unambiguously detect (at more 
than five statistical standard deviations) the presence of CP-violating phases over much of the 
allowed parameter space, although the effect diminishes with increasing \fj,\ and increasing tan /3 
(except in case B3, for the reason given above). For both these modes S and S give very similar 
results. Except for scenario Bl with strong Higgsino-gaugino mixing, pair production at 
both e^e~ and e~e~ colliders is much less promising, especially if the physical masses of Xi: Xi 
and Xg are held fixed, i.e. for <S. All these features can be understood from the discussion of 
Table 2 in Sec. 5.7. 

The small phase sensitivity of the e^e^ mode relative to the e^e^ mode can partly be 
explained by the smaller cross section of the former mode; recall that the significances scale 
with the square root of the number of events. In addition, closer inspection of the matrix 
elements shows that in case of e^e^ production, the terms oc cos0i and oc cos(0i + (p^) are 
suppressed by extra factors sin^ 9w and sin^ 9w relative to the leading phase-independent terms; 
for the e^e^ mode the corresponding relative factors are 1 and sin^ 9\y, respectively. 

Turning to chargino modes, we observe that they are sensitive to phases only in scenario 
B2, with large and for small tan/3. The only relevant phase here is 0^. Recall from the 
discussion of Sec. 4.3 that the maximal allowed value of this phase scales like This means 
that the maximal deviation of | cos0^| from unity scales like In case of xfxi production 
the main phase sensitivity comes from Tn^±, which gives an extra factor sin2/3/|//|. Altogether 
the maximal S{xi Xi ) therefore scales like |yLi|^ sin 2/3; this reproduces the numerical behavior 
in scenarios B2 and B3, with small gaugino-Higgsino mixing. A similar argument also holds 
for the mixed XiXt ^ode. However, in this case the cross section itself vanishes in the absence 
of gaugino-Higgsino mixing. This means that now the phase-dependent terms are of the same 
order in Mw as the phase-independent ones. Moreover, significant phase dependence now also 
comes from the Zxixt coupling, not only from the chargino masses. Hence both definitions 
of the significance now give very similar results. Finally, the very strong tan /3 dependence of 
these significances in case B2 is due to the fact that values of (p^ near tt are only allowed for 
small tan /3 in this case, see Fig. 2b. 

In contrast to the chargino modes, some neutralino modes are promising for all scenarios 
we considered. This is true in particular for the (12) mode. We saw in Sec. 5.7 that in this 
case both the total cross section and the phase dependence (on 0i) already start at 0{M^), i.e. 
they are not suppressed for large or large tan /3. Indeed, we find that this mode often allows 
somewhat better sensitivity than the celebrated e^e^ mode. The mixed gaugino-Higgsino 
modes also do well, especially for not too large values of |//|. As expected from the discussion 
of eqs.(5.33), the (Ih) mode is somewhat more promising than the {2h) mode. The rather good 
phase sensitivity of the (22) mode at first seems surprising, given that the phase dependence 
only enters at C(M|), whereas the cross section is 0{M^). However, closer inspection of the 
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sensitivities for the (22) and (Ih) modes shows that the relative factor between them is in fact 
0{\Mi\/Mz), which is close to unity in our case. Note that the relatively large size of the 
(22) cross section facilitates its precise measurements and therefore increases the significances. 
However, as remarked at the end of Sec. 6 we still consider the mixed (Ih) final state to be 
more promising, since it will be less sensitive to systematic uncertainties. 
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Table 4: The maximal significances S of eq.(6.2) and S of eq.(6.3) that can be found for choices 
of phases which are compatible with all low-energy constraints. The scenarios Bl, B2 and B3 
have been defined in Table 1. Notation and calculational procedures are as in Table 3. 



7.2 Polarizations 

As emphasized earlier, the significances S and S strictly speaking only measure deviations from 
the CP-conserving MSSM; they do not directly measure CP violation. Direct evidence for CP 
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violation could come from the measurement of the T-odd normal component of % polarization 
vectors introduced in Sec. 5.6. The maximal possible absolute values of these "polarization 
asymmetries" for scattering angle d — are summarized in Table 5. Recall that a nonzero 

asymmetry can emerge only in the production of two different x states, and that the asymmetry 
will be larger for the lighter of the two final-state particles. However, the polarization can only 
be measured through the % decay products; we therefore do not consider the polarization of 
X?, which is probably the LSP. 
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Table 5: Maximal absolute values of P^*^ in percent. The scattering angle 9 is set to |. 
Notations and conventions are as in Table 3. 

We see that the chargino polarization is likely too small to be useful, except in scenario 
B2 with large and small tan/3. Recall from the discussion at the end of Sec. 5.6 that this 
asymmetry (for the lighter chargino) scales like \n\ sin 2/3 sin0^; we saw in Sec. 4.3 that the 
upper bound on |sin0^| scales like Altogether the maximal value of P/v of the lighter 

chargino therefore scales like The very rapid decrease of this polarization with increasing 
tan P is partly due to the explicit sin 2/3 dependence, and partly due to the disappearance of 
the band around 0^ ~ vr, sec Fig. 2b. 

In scenarios with large \fi\ (B2, B3) the XiXl mode again proves most sensitive to CP- 
violating phases. Eq.(5.34) shows that in this case a nonzero Pn already emerges at 0{M^), 
and remains finite both for large and large tan/3. This describes well the behavior seen in 
cases where the perturbative diagonalization of the neutralino mass matrix is reliable. Moreover, 
recall from Table 3 that this mode has a fairly high cross section. This is important, since even 
for perfect (100%) analyzing power one needs nearly 1,000 events to detect a 10% asymmetry 
at the 3cr level. 

As expected from our earlier discussion of eqs.(5.37), the mixed gaugino-Higgsino modes 
also have sizable asymmetries even for large |//|, the heavier {2h) mode being more promising. 

However, the relatively small cross sections of these modes imply that one would need a very 
large luminosity for a meaningful measurement of polarization asymmetries in these modes. 
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except in scenario Bl with strong wino-Higgsino mixing. Indeed, in this last case the (13) and 
(24) modes are far more promising than the (12) mode. 

As noted earher, the spin of the produced x particles can only be determined on a statistical 
basis by (partly) reconstructing their decays. We find it encouraging that recent dedicated 
studies demonstrated sensitivity to phases in the ncutralino mass matrix using T-odd variables 
constructed in e~^e~ — >• XiXi with Xi Xi^^^^ [47], Xi ~^ f^r"^ — > t^t^Xi [48], and Xi ~^ Xi^ 
[49]. 

7.3 Correlations between observables 

In addition to their absolute sizes, the correlations between various phase-sensitive quantities 
are also of interest. Such correlations can provide stringent tests of the MSSM, since they are 
a consequence of the limited number of parameters affecting these leptonic observables in the 
MSSM. Recall that all our "high-energy" variables (cross sections and polarizations) depend 
on the phase 0^; most of them also depend on 0i, the exception being observables related 
to chargino pair production. We saw in Sec. 4.3 that 0^ is tightly constrained by the "low- 
energy" observables a^u and (especially) dg, while 0i in most scenarios can take any value (for 
some combination of the other phases). Moreover, the de constraint enforces a tight correlation 
between (J)^ and (pi, see Fig. 2. 

In Fig. 11 we compare high- and low-cncrgy quantities. We see that the phase-sensitive 
high-energy quantities are not correlated at all with d^. This is true both for T-even variables 
(Fig. 11a) and T-odd ones (b); in scenarios with strong gaugino-Higgsino mixing (a) and in 
scenarios where this mixing is suppressed (b); and for quantities that depend on both 0i and 
0^ (a) as well as those that depend only on 0^ (b). This can be explained from the observation 
made at the end of Sec. 4.3 that dg itself is not correlated with any of the phases after scanning 
over the other two phases; recall that the low-energy observables also depend on (pA- For 
example, except at the very edges of the allowed range of (p/^, de can still take any value within 
its experimentally allowed range even after is fixed; this is due to the variation of (pi and 

(pA- 

On the other hand, in some cases we do observe significant correlations between high- 
energy observables and a^. We saw in Fig. 5 that in scenarios B2 and B3 shows a cos— like 
dependence on (pi, in some cases (e.g. B2 at small tan/3) two separate bands of values exist, 
corresponding to cos (p^ ~ ±1. However, in scenario Bl shows very little correlation with (pi, 
see Fig. 5a, d. Correspondingly, Fig. 11c shows no correlation for scenario Bl, while Figs, lld-f 
show significant correlations for scenario B2. Comparison of panels d and e shows that this 
correlation becomes stronger at larger tan (3. This is due to the diminished role of (pA and the 
reduced width of the allowed band in the {(p^,(pi) plane; the overall size of |a^| also increases 
with increasing tan/3, see eq.(4.15). Finally, Fig. llf shows that high-energy quantities whose 
only phase sensitivity is through 0^ also correlate with a^. Note in particular that S{xtX2) 
is much bigger for < 0, which corresponds to 0^ ~ n, than for > 0, which corresponds 
to 10^1 ^ 1. This confirms the explanation we gave in the discussion of Table 4 for the very 
strong tan (3 dependence of this quantity. For this class of observables the correlation with a^u 
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Figure 11: Correlations between low- and liigti-energy quantities. Most high-energy observ- 
ables have been computed at y/s = 500 GeV, except for panels b) and f), which are for i/i = 800 
GeV. The parameter sets Bl and B2 have been defined in Table 1, and the significance S is 
defined via eq.(6.2). 



also becomes stronger with increasing tan/?; however, as remarked in Sec. 5.7, the sensitivity 
to 0^ disappears oc sin 2/5 at least. 

In most cases different phase sensitive high-energy observables are strongly correlated with 
each other. This is illustrated by Fig. 12, where we plot the two usually most promising 
significances, for the e^e^ and X1X2 final states, against each other. The simplest correlation 
obtains for scenario B3 for tan (3 = 20, shown in panel f. In this case the constraint excludes 
values of 0i near as well as 0^ near vr, see Fig. 2f. Hence the minimization in the definition 
(6.2) of S only goes over the single CPC point 0^ = 0, 0i = it. The strong correlation observed 
in Fig. 12f then follows from the fact that both significances shown here are essentially oc cos 0i 
to leading order in Mz, as explained in Sec. 5.7. 

The next simplest situation obtains if both 0i = and 0i = vr are allowed, but 0^ = vr is 
still forbidden, and tan/3 is not small [panels c), d) and e)]. Now the minimization in eq.(6.2) 
goes over two CPC points. Recall that this minimization is performed independently for the two 
significances shown in Fig. 12. The upper (lower) branch connected to the origin is populated 
by combinations of phases where both minimizations pick the CPC point 0i = (0i = vr). 
These two bands are connected by sets of points where our algorithm picks the CPC point 
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Figure 12: Correlations between tlie significances, defined as in eq.(6.2), for tlie processes 



e e 



e^Cj;^ and e'^e 



XiX2) botli taken at y/s 



500 GeV. 



01 = for 5(e^e^), but cfiooses tlie point 0i = vr for S{xiX2)- 

Fig. 12a shows that in scenario Bl the correlations get weaker at smaller tan/?. To un- 
derstand this, recall that scenario Bl has strong wino-Higgsino mixing, and hence a relatively 
strong dependence on 0^ through the combination cos(</)i + 0^), which depends linearly on 0^ 
when 1 01 1 and |0i — 7r| are sizable. In contrast, cos0^ depends only quadratically on 0^ for 
small 10^1, and can therefore to good approximation be set to 1 in scenario Bl, see Fig. 2a,d. 
This dependence on 0^ will be numerically different for the two modes present, loosening the 
correlation. This effect is important only at small tan/? for two reasons. First, all contributions 
to our cross sections that are sensitive to 0^ are suppressed by a factor sin 2/5 at large tan/3. 
Secondly, we saw that in scenario Bl the upper bound on |0^| decreases with tan/3. 

Fig. 12a shows another new effect on the lower branch, where both significances are evaluated 
with the CPC point 0i = tt, 0^ = 0. The cross section for X1X2 production in this case shows 
a non-monotonous dependence on cos0i. As expected from the expansion of the result (5.26) 
in powers of Mz using eqs.(2.17)-(2.19), this cross section reaches its absolute minimum at 
cos 01 = +1, where the S*— wave contribution vanishes. However, cos0i = —1 is also a (local) 
minimum, the maximum being reached at cos0i ~ —0.8; recall that the expansion in powers 
of Mz is not reliable in this case, since M2 = \fi\. As a result of this non-monotonous behavior, 
the cross section at cos0i ~ —0.6 becomes identical to that at cos0i = —1. Since (j(e^e^) 
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does decrease monotonically with cos 01, values of cos0i ^ —0.6 give rise to scenarios with 
very small S{xiX2) but sizable 5(e^e^). 

The comparison of Figs. 12b and e shows that the correlation becomes weaker for smaller 
tan/3 also in scenario B2. This is partly because the width of the allowed band in the {(j)^, </)i) 
plane decreases with increasing tan/3, see Fig. 2. In addition, in scenario B2 with tan/3 = 3 
the low-energy constraints also allow values of 0^ near tt. One can then find values of 0i not 
far from vr where cr{xiX2) CPV points with |0^| -C 1 is very close to this cross section 
at the CPC point ip^ = <pi = tt. This again leads to scenarios where S{x\X2) is very small, 
but iS(e^e^) is sizable. The existence of four different allowed CPC points also explains the 
occurrence of additional bands in Fig. 12b. 




10 20 30 40 50 

Figure 13: Correlation between the significances for the processes e~e^ ^l^l e~^e~ — 
Xixt, both taken at ^/s = 800 GeV, for scenario B2 with tan/? = 3. 

In some cases the correlations between different significances are quite weak. The most 
extreme case we found is shown in Fig. 13, and occurs for scenario B2 at tan/3 = 3. We saw in 
Table 4 that here (and only here) (t{xiX2) allows a significant probe of the phase 0^, whereas 
5(e^e^) is always mostly determined by 0i. Moreover, Fig. 2b shows that in the allowed band 
with 0^ ~ 71, the deviation |0^ — 7r| becomes maximal for 0i quite close to ±7r. This leads to 
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scenarios with large S^XiXt)^ but very small 5(e^e^). Conversely, || cos0i| — l| can be quite 
large for small \4>^\, leading to scenarios with 5(e^e^) ^ S{xixt)^ although the latter cannot 
be strictly zero if the former is bigger than 10. However, we saw earlier that other combinations 
of parameters do not allow meaningful probes of (f)^ using high-energy quantities. We therefore 
conclude that in most cases, significances that can be large are also fairly strongly correlated. 
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Figure 14: Correlation between the significance S and the absolute value of normal polarization 
Pat, measured at scattering angle 9 = n/2, for mixed neutralino pair production at ^/s = 500 
GeV. We consider X1X2 production for scenarios B2 and B3, but switch to the X1X3 state 
for scenario Bl. 



Finally, in Fig. 14 we compare the normal component of the polarization vector of the 
heavier neutralino in mixed neutralino pair production with the significance of the same mode. 
We consider X1X2 production in scenarios B2 and B3, but switch to X1X3 production in scenario 
Bl, where this final state is far more promising, see Tables 4 and 5. These figures look somewhat 
simpler than those in Fig. 12, since now the existence of two allowed CPC points only leads 
to two bands, as compared to three in Fig. 12. Of course, scenarios with a single allowed 
CPC point (Fig. 14f) again only yield a single band. In panel b) we again find scenarios with 

sizable phases, hence sizable \Pj^^' ^ ^|, and yet vanishing S{xiX2)'^ '^^ saw analogous behavior 
in Fig. 12b. 

More importantly. Fig. 14 shows that the polarization \Pn\ increases much more quickly as 
the (relevant) phase 0i is moved away from or vr than the significance S does. The reason is 
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that \Pn\, being T-odd, has a sin-hke dependence on 0i, i.e. grows hnearly with |0i| or — 7r|. 
In contrast, the T- and CP-even quantity <S has cos-hke dependence on all phases, and thus 
only grows oc or |0i — tt^ as 0i is moved away from a CPC point. T-odd observables like 
Pjv are therefore in principle better suited to probe small phases. 

8 Summary and Conclusions 

In this article we have discussed to what extent the phases of dimensionful parameters in the 
SUSY Lagrangian can be determined from leptonic observables. Since we assumed universal 
soft breaking parameters for the first two generations of sleptons and did not discuss processes 
involving third generation (s)particles, we only have to deal with three phases: those of the 
Higgsino mass parameter /x, of the f/(l)y gaugino mass Mi, and of the leptonic trilinear soft 
breaking parameter Ai, in all cases measured relative to M2 which we took to be real and 
positive by convention. 

Our main focus was on quantities that can be measured at future high-energy e'^e~ and 
e~e~ colliders, but we first analyzed the constraints that follow from the present measurements 
of the leptonic dipole moments de and a^. We worked in a scenario with moderately heavy 
sparticles; as well known, in this case sizable CP-odd phases are possible only if neutralino 
and chargino loop contributions to dg cancel to good approximation. In agreement with earlier 
work [10, 11], we found that, unless |//| » M2, m^-, the phases of Mi and Ai can take any value 
(for some combination of the other phases), whereas the phase of is tightly constrained, the 
maximal allowed deviation from or tt scaling like Our analysis of Sec. 4 also gave the 

perhaps surprising result that in this case improved measurements of d^ will not significantly 
reduce the allowed range for any one of the three relevant phases after scanning over the other 
two. This is true independently of whether this measurement leads to improved upper bounds 
on \de\ or finds a non-vanishing result. On the other hand, improved measurements of 
do have the potential to further restrict the allowed ranges of these phases; however, here 
improved measurements have to be combined with improved SM predictions for the hadronic 
contributions to a^. 

Turning to high-energy observables, we first analyzed in detail the phase sensitivity of total 
cross sections of various final states. To that end we introduced "significances" that determine 
the statistical significance with which the presence of non-trivial phases could be determined 
in a given production channel. As pointed out in ref. [38], the cross section for e^e^ production 
depends very strongly on the relative phase between Mi and M2; we found that a deviation of 
~ 60 to 90 standard deviations from the predictions of the CP-conserving MSSM is possible in 
this channel. However, this does not necessarily argue in favor of constructing an e~e~ collider, 
since certain neutralino production channels - in particular, X1X2 production for > M2 - 
have comparable or better sensitivity to the same phase. We also found a somewhat lower, but 
still promising, sensitivity in the e^ej final state. For our choice ~ 200 GeV, chargino pair 
production can show significant phase dependence over the experimentally allowed parameter 
space only for > 2M2. Since the de~constraint on 0^ becomes weaker for larger slepton 
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masses, the minimal ratio \\x\lM2 where chargino production channels can become useful for 
probing CP-violating phases should be smaller for larger m^-. However, these chargino modes 
will be useful only if tan/3 is quite small, since the relevant significances scale hke sin 2/3. 

A deviation of any of these cross sections from the prediction of the CP-conserving MSSM 
could perhaps also be explained by some extension of the model which does not introduce 
new CP-odd phases. We therefore also studied a CP-odd quantity: the component of the 
polarization of produced charginos and neutralinos that is normal to the production plane. We 
found that it can reach values exceeding 30% for the production of two different neutralinos; 
in scenarios with large \\x\ and small tan/3 the polarization vector of the lighter chargino, 
produced in association with the heavier one, could have an even larger normal component. 
Recent studies [47, 48, 49] indicate that such large CP-odd polarizations might indeed lead to 
measurable CP-odd asymmetries in the phase space distribution of the % decay products. 

Finally, we studied correlations between the various phase-sensitive observables. We found 
that the high-energy observables arc essentially not correlated at all with d^. This is due to 
the required rather precise cancellation between different contributions to d^] it implies that 
better measurements of d^ will not further restrict the possible ranges of phase-sensitive high- 
energy quantities. However, there is some correlation between these high-energy observables 
and a^. Moreover, most pairs of high-energy observables are quite strongly correlated with 
each other. This follows from the fact that most of them basically probe the phase of Mi, given 
the tight constraint on the phase of [i. Within the CP-violating MSSM the measurement of one 
phase sensitive high-energy observable therefore allows to greatly constrain the allowed range 
of other such quantities, thereby allowing stringent tests of the model. However, at large \\x\ 
and small tan /3 the phase of can play an important role, in particular in chargino production. 
In that case phase-sensitive observables in the chargino sector correlate poorly with those in 
the selectron or neutralino sector. This underscores the importance of measuring as many 
phase-sensitive quantities as possible. 

Total cross sections and CP-odd asymmetries offer complementary access to CP-odd phases, 
since they depend on these phases through cosine-like and sine-like functions, respectively. The 
former are rather insensitive to these phases if they are small (the perhaps most likely case). 
Measurements of, or bounds on, CP-odd asymmetries should then lead to better determina- 
tions or constraints on these phases. On the other hand, if some phase is near 7r/2, CP-odd 
asymmetries will be near maximal, which means that they are not well suited to precisely pin- 
ning down the value of this phase; precision measurements of some cross sections will then have 
the edge. Of course, there is also complementarity between high- and low-energy observables, 
since only the latter are sensitive to the phase of Ai. 

We conclude that measurements at high energy colliders will be necessary to pin down the 
phases of dimensionful parameters in the SUSY Lagrangian. Both precision measurements of 
CP-even quantities like masses and cross sections, and searches for CP-violating asymmetries, 
are promising in certain regions of parameter space. Linear e^e" colliders seem to be ideally 
suited for performing these measurements. 
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Appendix A Kinematics 

Working in the CMS frame with total energy a/s and neglecting the electron mass, the first 
electron and positron (second electron) momenta can be written as 

^5* = ^ (1,0, 0,1); (A.la) 

= ^ (1,0, 0,-1). (A.lb) 

The outgoing momenta of the produced superparticles b and c are: 

= ^ (l + ^^i^,Aisine,0,A|,cose^ , (A.2a) 

where A^c denotes the usual two-body final state kinematical function: 

A,e = A(^l,^,^); (A.3a) 

A(l, x,y) ^ 1 + x'^ + - 2{x + y + xy). (A.3b) 
Furthermore the kinematical invariants (Mandelstam variables) are 

s = (pi+P2)^; (A.4a) 

t = {Pi- kif] (A.4b) 

u = {pi- k2f. (A.4c) 



Appendix B Helicity amplitudes 

We calculate the relevant helicity amplitudes using the formalism introduced in [50] ^ Using our 
definition of tlie kiiieniatical situation, we find the following results for the scalar and vectorial 

^Our convention for a momentum-dependent Weyl spinor for fermions going in the —z direction differs by 
an overall sign from that of [50] . 
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fermionic string associated with massless fermions: 

Vip2, Cr2)PaUipi, (7i) = -aV^^aa/aio-a , (B.la) 
v{P2:Cr2h^Pau{pucri) = v^^^^j^^^ (0, 1, i(7i, 0), (B.lb) 

where the four choices in eq.(B.lb) correspond to n = 0, 1, 2, 3. In the case of neutrahnos or 
charginos with non-neghgible masses only the vectorial string is required. It can be written as 



Ui{h,Xi)^^Pf)Vj{k2,X2) = Y [Vl-^k'5AiA2(/9,Aisin^,0,AiCos^) (B.2) 

+ \J (1 + /?Air^^Ai)(l + /?Air^_;3Ai)(^Ai -A2(0,cos6', -iAi, - sin 9) 

where ^ 

Vpx^=Xl+f3Xl^^J, (B.3) 



and 



2 2 

m| — ruj 



A,, = (B.4) 



Appendix C Neutralino functions 

After introducing two effective neutralino mixing coefficients 

^ 2cos^w 2sin^w 

= (C.lb) 

^ cos^w 

we define two dimensionless neutralino functions for t- or u-channel exchanges: 

4 



M^^is, t/u) = J2 m^lV~sV^y^Dl^, (C.2a) 

k=l 
4 

N^p{s, t/u) = J2 ^VX^D^u, (C.2b) 



fe=i 



the propagators and D!^ have been defined in eq.(5.6). Very similar neutralino functions 
were introduced in [37]; we saw in Sec. 5 that they allow to give very compact expressions for 
the slepton production amplitudes. 
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